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A  paradigm  for  robust  geometric  algorithms 


1.  Introduction 

Geometrical  algorithms  in  use  in  computer  aided  design  systems  today  often  fail  for 
numerical  reasons.  The  cause  of  these  failures  usually  can  be  traced  to  logical  decisions  such 
as  'branch  on  zero"  that  depend  on  the  results  of  numerical  calculations.  Numerical  inaccura¬ 
cies,  introduced  either  in  the  initial  data  or  in  the  finite-precision  arithmetic  that  is  used,  may 
result  in  a  set  of  logical  decisions  that  are  inconsistent.  This  loss  of  logical  consistency  usually 
proves  fatal  to  the  algorithm 

Our  interest  is  in  geometric  algorithms  that  are  robust  in  the  sense  that  they  are  prov- 
ably  immune  to  such  potential  problems.  This  paper  explores  a  paradigm  that  should  have 
wide  applicability  for  producing  robust  algorithms,  and  it  applies  the  paradigm  to  the  task  of 
intersecting  two  convex  polyhedral  objects.  The  paradigm  was  previously  used  [HHK]  in  an 
implementation  of  an  algorithm  for  intersecting  polyhedral  objects  that  was  substantially 
more  robust  than  algorithms  implemented  earlier.  This  paper  is  a  first  step  in  developing  a 
mathematical  framework  that  justifies  the  underlying  ideas  in  this  implementation.  In  partic¬ 
ular,  an  important  tool  in  this  work  is  a  method  of  manipulating  e  mbedded  polyhedra  in  ways 
consistent  with  their  topology 

1.1  Correctness  paradigm 

The  computer  representation  of  a  geometrical  object  consists  of  two  types  of  data:  sym¬ 
bolic  and  numerical.  In  the  case  of  a  polyhedron,  the  symbolic  data  consists  of  names  of  ver¬ 
tices,  edges  and  faces  along  with  edge-vertex  and  face-edge  incidence  relations  Since  symbolic 
data  does  not  degrade  with  computation,  a  program  involving  only  symbolic  data  should  not 
fail  unless  there  is  a  flaw  in  the  program  itself. 

The  numerical  data  usually  consist  of  approximations  to  the  real  values  of  various  enti¬ 
ties.  In  the  case  of  a  polyhedron,  it  consists  of  the  coefficients  of  the  face  equations 


represented  to  some  precision.  For  a  program  that  performs  only  numerical  computations,  an 
argument  based  on  continuity  can  be  used  to  establish  that  the  result  is  an  approximation  to 
the  correct  answer.  Of  course,  such  a  program  may  fail  for  a  variety  of  obvious  reasons  such 
as  division  by  zero  or  an  attempt  to  find  the  square  root  of  a  negative  number  These  errors 
can  usually  be  handled  in  a  straightforward  manner.  A  more  subtle  type  of  error  occurs  when 
there  are  conversions  of  numerical  data  to  symbolic  data.  Typically  such  a  conversion  takes 
place  when  the  result  of  a  numerical  computation  is  used  in  a  branch  on  zero  operation.  If  the 
numerical  value  is  very  close  to  zero,  one  must  decide  if  the  quantity  is  actually  zero  with  the 
difference  due  to  round-off  error.  If  only  one  such  conversion  from  numerical  data  to  symbolic 
data  occurs,  a  slight  perturbation  of  the  numerical  input  data  should  force  the  result  of  the 
numerical  computation  to  agree  with  the  conversion.  Thus,  no  inconsistency  is  likely  to  arise 
from  either  decision.  However,  a  program  may  contain  a  sequence  of  such  conversions  that 
are  logically  dependent.  Thus,  inaccuracies  in  initial  data  or  in  numerical  computations  may 
result  in  inconsistent  conversions,  and  once  the  symbolic  data  is  inconsistent,  some  property 
which  is  essential  for  the  correctness  of  the  algorithm  may  fail  to  hold.  Indeed,  this  is  the  dif¬ 
ficulty  that  causes  many  implementations  of  geometrical  algorithms  to  fail.  The  solution  to 
this  problem  is  to  structure  the  algorithm,  so  that  all  conversions  are  logically  independent.  A 
major  contribution  of  this  paper  is  that  an  algc  '  .L  v  for  an  important  problem  such  as  the 
intersection  of  convex  polyhedra  can  be  structured  .ha  manner. 

1.2  A  simple  example 

Consider  the  following  algorithm  for  intersecting  two  polygons  Pi  and  P  2  Intersect 
each  edge  of  P{  with  each  edge  of  P Use  the  intersection  points  to  partition  edges  into  seg¬ 
ments.  Discard  each  edge  segment  of  P  \  that  does  not  lie  in  P  2  and  each  edge  segment  of  P 2 
that  does  not  lie  in  P , .  Assemble  the  remaining  segments  into  polygons. 

If  the  above  algorithm  is  applied  to  the  two  polygons  of  Figure  1,  numerical  round-off 
may  result  in  determining  that  edge  ( a,b )  intersects  edge  ( d,e )  in  some  point  h  but  that  edge 


(b,c)  fails  to  intersect  edge  ( d,e  i  In  assembling  edges  into  polygons,  it  is  likely  that  both  seg¬ 
ments  ( d,h )  and  (h,e)  have  been  discarded  as  lying  outside  P±.  In  this  case  there  is  no 
second  edge  at  vertex  h  to  associate  with  ( a,h ).  If  by  chance  (h,e)  is  classified  as  inside  P\, 
then  the  above  difficulty  will  surely  cause  problems  at  vertex  e.  In  either  case,  the  algorithm 
is  likely  to  return  a  structure  that  is  not  a  polygon 

In  this  simple  example,  it  is  clear  that  the  decision  as  to  whether  (a, 6)  intersects  (d,e) 
and  the  decision  as  to  whether  (b,c)  intersects  (d,e)  are  logically  dependent.  The  solution  to 
the  difficulty  is  to  modify  the  computation  so  that  all  conversions  from  numerical  data  to  logi¬ 
cal  data  are  logically  independent.  In  the  above  example,  this  can  be  done  by  determining  on 
which  side  of  (d,e)  b  lies,  and  then  using  symbolic  reasoning  to  determine  whether  ( a,b )  and 
(d,e)  intersect  and  whether  (6,c)  and  (d,e)  intersect. 

The  range  of  problems  for  which  algorithms  can  be  written  so  that  all  conversions  from 
numerical  data  to  symbolic  data  are  independent  is  quite  broad  and  as  we  shall  see  includes 
the  intersection  of  convex  polyhedra.  It  is  important  to  determine  the  extent  of  the  class  of 
problems  for  which  this  can  be  done.  Problems  outside  the  class  can  be  classified  by  the  power 
of  a  theorem  prover  needed  to  establish  that  a  given  set  of  conversions  is  or  is  not  indepen¬ 
dent 


Figure  1  Intersecting  two  polygons. 

1.3  Definition  of  correctness 

A  definition  of  correctness  for  geometric  algorithms  must  take  into  account  the  fact  that 
numerical  input  data  only  approximates  the  real  date  and  thus  may  be  inconsistent  with  the 
logical  input  data.  The  following  definition  of  correctness,  adopted  in  [HHK],  allows  for  this 
possible  input  inconsistency.  It  depends  on  suitable  notions  of  approximation,  computer 
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representation,  and  geometric  entity  or  object,  which  we  discuss  in  §14. 

We  say  that  an  algorithm  for  an  operation  °  is  correct  if,  for  any  representations  R  i , 
R  2.  '  '  it  produces  a  representation  R  such  that  there  exists  entities  M  M 2,  and  M 

such  that 

(1)  the  numerical  data  given  in  R  and  Ri  are  approximations  to  the  numerical  data 

given  in  M  and  M, , 

(2)  the  logical  data  in  R  and  Rt  agree  with  those  in  Af  and  Mlt 

(3)  2.  ■  ■  ■  )  =  M. 

Note  that  correctness  of  an  algorithm  in  the  above  sense  is  not  sufficient  when  the  algo¬ 
rithm  is  used  as  a  subroutine  in  a  more  general  algorithm.  Consider  an  algorithm  for  inter¬ 
secting  line  segments  that  reports  an  intersection  only  when  the  line  segments  themselves  or 
a  slight  perturbation  intersect,  and  reports  a  nonintersection  only  when  the  line  segments 
themselves  or  a  slight  perturbation  do  not  intersect.  Such  an  algorithm  is  sufficient  for  a  sin¬ 
gle  instance  of  testing  whether  two  line  segments  intersect.  However,  the  same  algorithm 
may  not  be  satisfactory  in  a  more  general  setting.  For  example,  the  edge  intersection  algo¬ 
rithm  in  the  earlier  polygon  intersection  algorithm  may  be  correct  by  the  above  definition  but 
fail  to  work  satisfactorily  in  the  polygon  intersection  setting. 

1.4  Computer  representations 

A  geometrical  object  such  as  a  polyhedron  is  regarded  as  a  mathematical  structure  con¬ 
sisting  of  vertex-edge  and  edge-face  incidences  along  with  equations  for  the  faces.  The  equa¬ 
tions  for  the  faces  have  coefficients  that  are  real  numbers.  The  geometric  incidences  deter¬ 
mined  by  these  equations  match  the  listed  vertex-edge  and  edge-face  incidences.  A  computer 
representation  or  3imply  a  representation  of  the  object  is  a  data  structure  presenting  the 
incidence  relations  and  face  equations.  However,  in  the  computer  representation  coefficients 
of  the  face  equations  are  approximated  by  finite  precision  numbers.  The  vertex-edge  and 
edge-face  incidences  implied  by  the  given  face  equations  need  not  agree  completely  with  the 
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incidence  information  in  the  representation  In  particular,  where  the  incidence  information 
requires  four  or  more  planes  to  meet  at  a  vertex,  the  actual  geometrical  data  may  imply  that 
the  planes  meet  in  a  more  complex  structure  consisting  of  several  vertices  connected  by  a  tree 
of  short  edges  icf.  §4.1)  The  incidence  information  requires  this  structure  to  collapse  to  a  sin¬ 
gle  vertex,  thus  imposing  a  condition  on  the  faces  that  is  not  stable  under  small  perturbation 
of  the  face  equations.  Representations  are  discussed  further  in  §4.1. 

1.5  Conceptual  overview 

In  this  paper  the  paradigm  for  robust  computation  is  applied  to  the  problem  of  intersect¬ 
ing  two  convex  polyhedra.  The  intersection  is  carried  out  by  intersecting  one  of  the  polyhedra 
with  each  half  space  defining  the  other.  A  plane  defining  a  half  space  may  cut  the  first 
polyhedron  so  that  certain  vertices  are  on  the  plane  or  are  so  close  to  the  plane  that  we  cannot 
numerically  decide  whether  they  are  or  are  not  on  the  plane.  In  this  case  we  will  declare  each 
such  vertex  to  lie  on  the  plane.  We  need  to  prove  that  we  can  slightly  perturb  the  face  equa¬ 
tions  of  the  polyhedron  so  that  the  vertices  claimed  to  be  on  the  plane  indeed  are  on  the  plane. 
The  difficulty  here  is  that  faces  incident  to  a  high-degree  vertex  may  no  longer  meet  in  a  point 
after  the  perturbation.  Thus,  we  need  to  develop  tools  that  will  allow  us  to  perturb  face  equa¬ 
tions  of  polyhedra  without  creating  small  features.  One  of  these  tools  is  an  isomorphism 
between  certain  stressed  planar  graphs  and  convex  polyhedra. 

Consider  a  planar  graph  whose  exterior  vertices  are  fixed  to  form  a  convex  polygon  and 
whose  interior  edges  behave  like  springs.  The  interior  vertices  will  come  to  rest  at  an  equili¬ 
brium  position  where  the  forces  at  each  interior  vertex  sum  to  zero.  Suppose  that  we  can  then 
apply  stresses  to  the  edges  of  the  exterior  polygon  so  that  the  forces  at  the  exterior  vertices 
also  sum  to  zero  The  graph  is  then  said  to  be  in  equilibrium.  There  is  a  correspondence 
between  such  graphs  in  equilibrium  and  convex  polyhedra.  It  is  easier  to  perturb  the  vertices 
of  the  graph  by  adjusting  the  spring  constants  and  thereby  perturbing  the  location  of  the  ver¬ 
tices  of  the  associated  polyhedron  than  it  is  to  perturb  the  vertices  of  the  polyhedron  directly. 
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In  Section  3  we  develop  the  correspondence 

Certain  algorithms  that  deal  with  polyhedra  require  that  the  polyhedron  have  no  small 
features;  e  g.,  two  vertices  extremely  close  together  Such  is  the  case  with  our  algorithm  for 
intersecting  a  polyhedron  with  a  half  space.  However,  intersecting  a  polyhedron  with  a  half 
space  may  produce  a  polyhedron  with  small  features.  Thus,  to  iterate  the  intersection  algo¬ 
rithm  we  need  a  method  to  remove  small  features  after  each  intersection  operation.  Two 
options  for  removing  small  features  are  discussed  in  §4.5 . 

In  Section  5  we  apply  the  isomorphism  of  Section  3  to  the  problem  of  intersecting  a  con¬ 
vex  polyhedron  with  a  half  space.  For  technical  reasons,  the  construction  and  proof  are 
simpler  in  the  special  case  where  the  polyhedron  has  a  triangular  face  and  the  plane  defining 
the  half  space  cuts  this  triangular  face  into  two  reasonable  size  pieces.  The  general  case  is 
briefly  described  in  5.1. 

1.0  Statements  of  results 

There  are  four  major  contributions  in  this  paper. 

•A  paradigm  for  constructing  robust  geometrical  algorithms. 

•A  demonstration  that  decisions  based  on  numerical  calculations  can  be  made  in  a 
consistent  manner  fcr  a  nontrivial  problem:  namely  the  intersection  of  a  convex 
polyhedron  and  a  half-space. 

•An  isomorphism  between  weighted  graphs  and  polyhedra  that  can  be  used  to  mani¬ 
pulate  the  geometrical  embedding  of  polyhedra  in  ways  consistent  with  their  topol¬ 
ogy 

•A  robust  algorithm  for  intersecting  a  convex  polyhedron  and  a  half  space  that 
should  never  fail  for  numerical  reasons  Under  suitable  assumptions  on  precision, 
the  algorithm  iB  correct  in  the  sense  of  §1.3. 
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2  Background  on  convex  stressed  graphs 

In  this  section  we  introduce  the  concept  of  an  equilibrium  stress  on  a  planar  graph  fol¬ 
lowing  definitions  in  [C]  and  [W],  Given  an  abstract  planar  graph  with  a  set  of  positive 
stresses,  together  with  a  peripheral  cycle  of  the  graph  embedded  as  a  convex  planar  polygon, 
there  is  a  unique  extension  to  a  convex  embedding  of  the  entire  graph  satisfying  a  given 
equilibrium  condition  at  each  interior  vertex  (cf.  Theorem  2.1.3  and  Theorem  2.1.4).  The 
importance  of  this  result  is  that  there  is  a  one  to  one  correspondence  between  sets  of  positive 
stresses  and  planar,  convex  embeddings  of  the  given  graph  with  fixed  periphery. 


2.1  Equilibium  stresses  on  planar  graphs 

Let  T  denote  an  abstract  planar,  3-connected  graph.  The  set  of  such  graphs  coincides 
with  the  set  of  1-skeleta  of  convex  polyhedra  [G,  p. 235].  A  mapping  of  the  n  vertices  of  T  to 
points  in  the  plane  H  =  {(x.y.z)  |  z  =  11  C  R  3  determines  a  realization  of  T  Let  us  denote 
the  n-tuple  of  these  points  p  =  (pi,p2.  '  '  and  the  realization  T(p).  If  the  points  are 

distinct  and  the  open  edges  disjoint,  then  we  say  that  T(p)  is  embedded.  An  embedded  graph 
T(p)  is  convex  if  its  faces  are  convex  and  its  periphery  iB  a  convex  &-gon.  Edges  not  on  the 
periphery  are  called  interior  edges.  It  will  be  useful  to  select  a  fixed,  abstract,  peripheral  k- 
gon  in  T  (  cf.  [  T]  i  and  to  denote  it  by  11.  In  realizations  Hp)  of  T,  we  may  denote  the 
corresponding  realization  of  IT  by  Il(p).  We  assume,  unless  stated  otherwise,  that  vertices  are 
numbered  so  that  P\,P2-  '  '  '  >Pk  correspond  to  the  vertices  of  ri(p). 

A  stress  on  T  is  a  collection  («,,}  of  real  numbers,  i  ,j  =  1, 2,  •  •  ,  n ,  t  *  J,  satisfying 

a  symmetry  and  a  vanishing  condition: 


(a)  uly-wyi,  for  all  i,j,  t  * j 

ib)  u,j  =0,  unless  p,  is  connected  to  Pj  by  an  edge. 


(1) 


If  each  to ,  corresponding  to  an  edge  of  T  is  strictly  greater  than  zero  or  less  than  zero  accord¬ 
ing  as  the  edge  is  interior  or  peripheral,  we  say  that  is  convex. 
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For  a  given  choice  of  p  -  {P\P  x-  '  '  '  .p,).  the  so-calied  equilibrium  condition  at  p. 

i«yip.-p,->  =  ° 

;  =1 

may  be  satisfied  for  all  pL,  in  which  case  we  call  {to,,}  an  equilibrium  stress  and  '{utJ\,  V(p ))  a 
stressed  graph.  If  i2>  is  satisfied  only  when  p ,  is  an  "internal”  vertex,  we  call  i <*),_,}  a  restricted 
equilibrium  stress  and  ({w,,},  T(p))  a  restricted  stressed  graph.  For  convenience,  we  also 

define  a  restricted  convex  stress  on  T  to  be  a  stress  {<•),;}  on  T  such  that  each  internal  w,,  is 
strictly  greater  than  zero  with  no  condition  on  the  peripheral  stresses. 

Recall  that  the  vertices  p,  belong  to  H,  so  that  in  standard  coordinates  they  are  of  the 
form  p,  =  (x,,_yt,  1).  Sometimes  it  will  be  convenient  to  use  instead  the  coordinates  (x..^,). 
This  clearly  results  in  no  essential  change  of  information. 

Following  Connelly  [C],  we  start  with  a  fixed  restricted  convex  stress  for  T,  and  we 
define  an  associated  energy  function. 

E(p)  =  )r  2  UlJ\pl-pJ\2,  (3) 

where  p  =  (pi.p2.  ,p„)  ranges  over  Hn.  Clearly,  E  is  a  homogeneous  quadratic  func¬ 
tion  of  the  coordinates  x  t ,  y  i ,  •  ■  ,  x„ ,  yn . 

We  now  want  to  hold  Xi.yi,  •  •  ,x*,y*,  (  i.  e.,  p i,  ‘  p*)  fixed.  Thus,  E  becomes  a 

nonhomogeneous  quadratic  function  of  the  remaining  variables:  E(p )  —  +  .yn) 

Set  e  =  u*  +i,  '  ’  •  ,y„).  Then  we  may  write 

F (z)  =  Qiz)  +  L(z)  +  Fi 0)  , 

where  Q(z)  (resp.,  Liz),  F{ 0))  is  the  sum  of  the  quadratic  (rasp.,  linear,  constant;  terms  of  F 

Define  a  bilinear  form  Biz  ,z  )  by  the  equation 

Biz  .z  )  =  Qiz  +  z  )  -  Qiz  )  ~  Qiz  ). 

Set  m  =  n  -  k  and  define  the  2m  X  2m  matrix  B  b>  the  rule 


B[z,w)  =  u>‘  B  •  z. 
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Similarly,  <*-■'  ie  the  ‘column1  '2m  -  tupie  L  by 

L  ( u' )  =  iv‘  L. 

Fhen  the  gradient  of  F  is  given  by 

VF(z)  =  B  ■  z  +  L. 

From  equation  (3>  above,  we  can  immediately  compute  that  the  l2i  —  1  >>r  and  2i!h  coor- 
dinates  of  + 1 ,  ,p„)  are  given  together  by 

n 

2  u,j(p,  -  pj),  for  i  =  k  +  1,  •  •  •  ,n. 
j  =i 

Thus,  we  may  immediately  conclude  the  following: 

2.1.1  Lemma:  The  following  are  equivalent: 

la)  | u.j)  is  a  restricted  equilibrium  stress  for  T(p). 

ib)  z  =  (pk  +l ,  •  •  ,pn)  is  a  critical  point  of  F  {  i.  e.,  VF(p*  + 1 ,  ,p„)  =  0). 

(c)  B  ■  z  =  -L.  □ 


We  note  here  that  after  a  suitable  permutation  of  rows  and  columns,  the  2mX2m 
matrix  B  actually  has  the  simpler  form 

C  0 
0  C 

for  a  certain  mXm  matrix  C  (see  Lemma  2.2.1(b)) 

2.1.2  Lemma.  F  has  a  unique  critical  point,  at  which  it  achieves  its  minimum 

Proof:  The  proof  is  modeled  on  a  proof  in  Connelly  [C].  Suppose  that  z  =  (p*  +  ;,  ,pn)  is 
large  Then  p,  is  large  for  some  j  >k.  This  impues  that  in  the  shortest  edge-path  connecting 
Pj  to  some  p, ,  i  S  k  ,  some  edge  is  large.  But  the  coefficient  in  (3)  corresponding  to  this  edge 
is  strictly  greater  than  zero,  and  so  £(p)  =  F{z)  is  large  Thus,  for  sufficiently  large 
a  >  0,  |  >a  implies  that  F(z)  >  F( 0)  It  follows  that  F  must  achieve  its  minimum  on 
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(2  I  |*  |  ^a} 

To  prove  uniqueness,  suppose  F  achieves  its  minimum  at  z'  and  that  z'1  is  a  critical 
point  of  F  For  real  t  ,  define  fU)  =  F{tc'  +  (1  —  t)z").  This  function  is  at  most  quadratic 
in  t.  Hence,  its  derivative  f'(t)  is  at  most  linear.  By  assumption.  f'(C)  satisfies 
f'( 0)  =  f'{  1)  =  0.  Thus,  f  is  constant  Bui  then  tz'  +  (1  —  t)z"  must  be  constant  since 
F(z)  is  large  for  large  z.  Therefore  z"  =  z' .  Q 

2.1.3  Theorem:  Choose  a  convex  k-goa  in  H  with  consecutive  vertices  q  ! ,  qr  2  •  '  .  <?* •  Sup¬ 

pose  {ulj}  is  any  restricted  convex  stress  on  T.  Then  there  exists  a  unique  realization  Hp)  of 
T  satisfying: 

(1)  p,  =<?,,<  =  1,2,  •  •  •  ,k, 

(21  {<*>,,}  is  a  restricted  equilibrium  stress  for  Tip),  and 

(3)  the  points  Pk+\.Pk+2<  '  <Pn  depend  smoothly  on  the  ui;,  i  +  1  or  j  3:k  +  1. 

Proof:  Define  the  energy  function  E  as  in  (3).  Set  p,  =  q,,i  =  1,2,  ,k,  and  define  F 

as  above.  By  Lemma  2.1.2,  F  has  a  unique  critical  point  z  =  (p*+i,  •  ,p„),  which 
implies,  by  Lemma  2.1.1,  that  is  a  restricted  equilibrium  stress  for  Hp),  as  desired. 

Uniqueness  of  the  realization  follows  from  the  uniqueness  portion  of  Lemma  2.1.2. 

Finally,  to  see  that  z  =  (,p*+l,  •  ,p„)  depends  smoothly  on  |u, ,  1 1  S  k  +  1  or 

j  £  k  +  1},  we  consult  equation  (c)  of  Lemma  2.1.1:  Bz  =  —  L.  As  we  shall  see  in  Lemma 
2.2.1  below,  B  and  L  depend  smoothly  on  the  above  stress  constants.  To  complete  the  proof, 
then,  it  remains  to  show  that  the  matrix  B  is  invertible,  for  then, 

z  =  -(B  _l )  L, 

and  smooth  dependence  of  z  follows. 

But  the  invertibility  of  B  follows  immediately  from  the  equivalence  2.1.1(b)  <=>  2.1.1(c), 
together  with  the  uniqueness  part  of  Lemma  2.1.2.  Q 


Tutte  [Tj  obtains  a  special  case  of  Theorem  2  1.3,  making  use  of  an  inductive  computa¬ 
tion  in  an  earlier  paper  [BSST],  The  proof  of  the  general  case  presented  here  is  based  on  some 
ideas  of  Connelly  [C]  and  leads  nicely  to  subsequent  calculations. 

2.1.4  Theorem  (  Tutte  [T]  ):  Suppose  that  T(p)  is  a  realization  of  T  such  that  [l(p)  is  a  con¬ 
vex  k-goa.  Suppose  that  there  exists  a  restricted  convex  stress  on  T  which  is  a  restricted 
equilibrium  stress  for  Tip).  Then  T(p)  is  convex.  Q 

In  particular,  the  existence  of  a  conu°x,  restricted  equilibrium  stress  for  T(p)  <  with  F  3- 
connected)  implies  that  T(p)  is  embedded,  with  all  nonperipheral  vertices  pk  +  i ,  •  ,p„  actu¬ 

ally  interior  to  ri(p).  Tutte  [T]  proves  this  result  in  the  case  that  all  the  interior  equal 
one,  but  the  proof  extends  virtually  verbatim  to  our  more  general  case. 

Theorem  213  tells  us  that  if  we  pin  the  peripheral  vertices  at  pi,p2,  ,p*  so  as  to 
form  a  convex  k-gon  in  H,  then  any  assignment  of  positive  stress  constants  to  internal  edges 
of  T  will  uniquely  and  smoothly  produce  internal  vertices,  at  each  of  which  the  equilibrium 
conditions  are  satisfied.  Theorem  2.1.4  asserts  that  the  resulting  T(p)  is  convex  (hence  embed¬ 
ded)  Unfortunately,  Theorem  2.1.3  (2)  above  gives  us  only  a  restricted  equilibrium  stress, 
whereas  our  later  applications  require,  in  addition,  equilibrium  at  each  peripheral  vertex 
Pi  P2,  <Pk-  Theorem  2.3.1  below  shows  that  we  can  obtain  these  additional  equilibrium 
conditions  when  k  =3. 


2.2  Some  computations 

Given  a  stress  {«,_,}  on  I\  we  follow  [C  ]  and  define  an  associated  stress  matrix  Q  as  fol¬ 
lows: 


0, 


—  w 


u  • 


<  *J 


n 


2 

*  =  i 


*  =j 


t,j  =  1.2.  •  •  •  ,n 


i.4 1 


The  equilibrium  equations 
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2  w,,tp,  -  Pj)  =  0,  i  =  1,2,  ■  •  •  .  n 

j  =i 

may  now  be  rewritten  as  the  pair  of  matrix  equations 

G  •  X  =  0 
G  •  Y  =  0 

where  X  (resp,  Y  )  is  the  column  n-tuple  (x!.X2,  •  ,xn)‘  (resp.,  (yi,V2. 

ing  of  the  x-coordinates  (resp.,  y-coordinates)  of  the  p,  =  U,,y,.  1). 


(5) 


(6) 

,y„)‘)  consist- 


Recall  that  vertices  have  been  numbered  so  that  the  vertices  of  the  peripheral  k-gon 
nip)  are  pi,P2.  '  '  '  P*  Then,  we  write  G  in  block  form  as  follows 


G  = 


A  B 
B‘  C 


(7) 


where  A  is  k  Xk  and  C  is  m  X  m ,  m  —  n  -k.  It  is  easy  to  see  that  the  equilibrium  equa 
tions  (5)  hold  for  t  =  k  +  1,  k  +  2,  ■  •  ■  ,n,  if  and  only  if 


B‘X  +  CX  =  0 
B‘Y  +  CY  =  0 


(8; 


where  X'  =  (x  i ,  ■  ■  ■  ,xk)‘ .  X  =  (x*  +  l,  •  •  •  ,x„)‘  and  similarly  for  Y  and  Y  Thus,  (8) 
expresses  the  fact  that  is  a  restricted  equilibrium  stress  on  T(p). 

Comparing  i&>  with  Lemma  2.1.1,  we  can  easily  conclude  the  following: 


Lemma  2.2.1:  <a>  Up  to  a  permutation  of  coordinates, 


B‘X' 
B‘Y' 

(b)  Up  to  a  permutation  of  rows  and  columns 


B  = 


C  0 
0  C 


□ 


Corollary  2.2.2:  Assume  that  pl}  ■  ■  ,Pk  (hence  A '  ,Y  )  are  fixed,  and  that  is  a  res¬ 
tricted  convex  equilibrium  stress  on  T.  Then  C  is  invertible,  and  the  vertices  p*  +  i,  '  .p* 
produced  by  Theorem  2.1.3  (  i.  e.,  X  ,Y  )  may  be  computed  by 
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X  =  —C  ~lB‘X' 

(9i 

r  =  -C~lB‘Y 

Proof:  The  proof  of  Theorem  2.1.3  shows  that  B  is  invertible.  Hence  C  is,  by  Lemma  2.2  1(b) 
above.  Then  <9)  follows  immediately  from  (8).  Q 


2.3  Finding  stress  constants  for  graphs  with  triangular  periphery 

2.3.1  Theorem:  Suppose  that  Tip)  is  a  realization  of  T  with  Tl(p)  a  nondegenerate  triangle 
(i.  e  ,  k  =  3).  Suppose  also  that  {<jiy}  is  a  restricted  equilibrium  stress  for  Hp).  Then,  there 
exists  a  unique  equilibrium  stress  for  T(p)  satisfying: 

1 1)  u,y '  =  ulJt  if  i  or  j  >  3. 

(2)  (<j,y'|l  2S  i,j  S  3}  depends  smoothly  on  («ly  1 1  or  j  >  3}. 

(3)  If  |w,y}  is  restricted  convex,  then  {<oly’}  is  convex. 


Proof:  We  use  the  notation  and  computations  in  2.2.  Thus,  {uly}  corresponds  to  a  stress 
matrix 


fl  = 


A  B 
B‘  C 


as  in  (7),  and  the  x-  and  y-coordinates  of  p  yield  column  n-tuples 

X  =  ]  and  y  =  [  y  ',' ) 

where  X'  and  Y'  are  column  triples,  etc.  Let  1  denote  the  column  n-tuple  (1,1,..., I)1,  with 
I'  and  1"  the  corresponding  triple  and  (n  —  3)-tuple. 

Let  Q;  denote  the  stress  matrix  corresponding  to  the  sought  after  equilibrium  stress 
{<d,y Condition  (1)  of  the  theorem  requires  that  Q'  has  the  form 

f  A'  B  1 


where  A'  is  a  3X3  matrix  to  be  determined.  The  equilibrium  condition  is  equivalent  to  the 
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equations 

G'  A'  =  0 
Q'Y  =  0 

whereas  the  fact  that  Q'  has  row-sums  equal  to  zero  is  equivalent  to 

o'i  =  o. 

These  yield  the  following  equation  involving  A 

A'[X\  Y',V]  +  B  [X",Y", I"]  =  0. 

Since  the  triangular  periphery  of  Hp)  is  nondegenerate,  the  3X3  matrix  [A'',  is  non¬ 

singular.  Thus,  the  last  equation  determines  A  '  uniquely: 

A '  =  -BiX",Y",V,HX',Y',V]-1. 

This  demonstrates  the  existence  and  uniqueness  of  an  equilibrium  stress  {<*>,/}  on  T(p) 
satisfying  (1)  and  (2).  It  remains  to  demonstrate  (3),  which  we  do  by  a  geometric  argument. 

Note  that  if  equilibrium  is  satisfied  at  a  vertex,  there  cannot  exist  a  line  through  that 
vertex  separating  the  positively  stressed  edges  from  the  negatively  stressed  ones.  It  follows 
immediately  that,  under  conditions  of  equilibrium,  if  the  internal  edges  meeting  a  given  peri¬ 
pheral  vertex  are  all  positively  stressed,  then  the  peripheral  edges  meeting  that  vertex  must 
be  negatively  stressed.  Q 

3.  Polyhedra  and  stressed  planar  graphs 

In  this  section  we  describe  a  basic  correspondence  between  certain  planar  graphs  and  cer¬ 
tain  polyhedral  surfaces.  This  correspondence  allows  us,  in  some  cases  of  interest,  to  reduce 
problems  of  manipulating  polyhedra  to  problems  of  manipulating  graphs.  Our  view  is  that 
this  gives  us  a  modeling  tool  with  potentially  wide  applicability.  We  back  up  this  claim  in 
Section  4  with  applications  to  removing  small  features  and  intersecting  polyhedra  with  half 


spaces. 
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The  families  of  polyhedral  surfaces  and  of  stressed  planar  graphs  that  we  consider  have 
natural  structures  as  vector  spaces,  with  one  mapped  to  the  other  by  a  linear  isomorphism 
The  mapping  was  first  noted  by  Maxwell  [M]  and  the  bijectivity  property  by  Crapo  and  White- 
ley  [W],  Our  description  of  the  mapping  differs  somewhat  from  earlier  descriptions  and  has 
several  advantages  It  is  given  by  an  explicit  formula,  it  is  visibly  linear,  and  it  uses  only 
well-known,  elementary  linear  algebra  The  proof  of  bijectivity,  however,  requires  a  homologi¬ 
cal  argument.  This  follows  easily  from  our  formulation  and  is  given  in  Appendix  B.  In  fact, 
this  argument  is  a  special  case  of  a  quite  general  homological  construction  that  will  appear 
elsewhere  [K]. 

The  authors  wish  to  thank  Robert  Connelly  for  a  number  of  helpful  conversations  and 
suggestions  concerning  the  material  in  this  section. 

3.1  Stressed  graphs 

Let  Hp)  be  a  3-connected,  planar  graph  embedded  in  the  plane  H  =  {(x,y,*) |z  =  1}, 
with  boundary  k-gon  II(p).  Let  S(I\p)  denote  the  set  of  equilibrium  stresses  (to,,)  on  T(p). 
We  may  define  addition  and  scalar  multiplication  of  equilibrium  stresses  by  the  rules 

{«„}  +{«.;}  =  + 

c 

With  these  operations,  SlT.p)  becomes  a  vector  space. 

3.2  Polyhedral  surfaces 

Let  fi,  f%,  •  •  •  denote  the  interior  faces  of  Tip),  and  let  /©  be  the  closed  region  bounded 
by  nip),  which  we  also  call  a  face.  We  now  describe  polyhedral  surfaces  that  "sit  over"  the 
graph  Tip).  Intuitively,  such  a  surface  may  be  constructed  from  planar  polygonal  faces  in  R3, 
one  for  each  face  in  Tip)  Each  face  in  R3  projects  orthogonally  down  to  its  counterpart  in 
T(p)  and  matches  up  on  edges  with  neighboring  faces.  In  order  to  deal  with  these  surfaces 
algebraically,  we  recognize  that  if  fr  is  a  face  of  Tip),  then  any  polygonal  face  in  R3  that  pro- 
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jects  to  fr  may  be  envisioned  as  the  graph  of  an  affine  function  defined  on  fr.  Thus,  a 
polyhedral  surface  may  be  represented  as  a  certain  family  of  affine  functions.  We  now  make 
this  precise 

Recall  that  Tip)  sits  in  the  plane  H.  Let  A(H)  denote  the  set  of  all  affine  functions 
H—*R  A  typical  such  function  a  is  given  by  a  formula  a(x,v,  l)  =  cx  +  d\+e,  for  certain 
constants  c,d,e.  Addition  and  scalar  multiplication  of  functions  make  A(H)  into  a  vector 
space  of  dimension  three.  Let  jar}  be  a  sequence  of  affine  functions  in  A  (H),  one  for  each  face 
fr  of  T(p)  (including  /'o).  We  say  that  the  sequence  {ar}  is  ptecewise-affinr  (P.A.)  if  the  ar 
satisfy  a  suitable  compatibility  property  Namely,  for  each  edge  etJ  and  incident  faces  fr,fs, 
we  must  have 

ar  I  etJ  =as  |  e,j  (10) 

where  a  |e  is  the  restriction  of  a  to  e.  If  we  imagine  the  graphs  of  ar  \  fr  and  as  |  fs  as  sitting 
in  R3  "over"  the  plane  H,  then  <L0)  simply  announces  that  these  graphs  match  up  over  the 
edge  etJ.  Thus,  the  sequence  {ar\  is  a  representation  of  a  surface  sitting  over  f(p). 


Figure  2:  Polyhedron  sitting  over  plane  and  the  graph  of  its  projection. 

We  now  collect  all  the  piecewise-affine  sequences  { ar }  into  a  set,  denoted  PA(Y,p),  and 
note  that  the  usual  operations  of  termwise  addition  and  scalar  multiplication  of  sequences 
make  PA  (T.p)  into  a  vector  space.  By  our  remarks  above,  this  vector  space  may  be  identified 
with  the  space  of  polyhedral  surfaces  over  Hp). 

Now  suppose  that  {6r}  is  a  P  A.  sequence  and  6o  =  0.  This  means  that  in  the  correspond¬ 
ing  polyhedral  surface,  the  face  represented  by  the  graph  of  b0\f0  actually  coincides  with  f0. 
In  effect,  the  polyhedral  surface  "sits  on"  fo  (in  H)  If  {ar}  is  any  P.A.  sequence,  we  can  sim¬ 
ply  subtract  qq  from  each  ar  to  obtain  a  P.A.  sequence  sitting  on  [q.  This  translates  into  the 
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following  simple  assertion  about  the  vector  space  PA(T.p)  of  all  P  A  sequences.  Let 
PA(r,p)0  consist  of  all  P.A.  sequences  {ar}  sitting  on  < i.e . ,  a q=0>.  This  is  a  vector  sub¬ 
space  of  PAlT.p),  as  is  the  set  of  all  constant  sequences  (i.e.,  sequences  {ar}  with  ar=ao,  for 
all  r).  Together,  these  subspaces  give  a  direct  sum  decomposition  of  PA  (V ,p): 

PA(V.p)  =  PA(T,p) o  +  /constant  sequences/. 

Since  a  constant  sequence  is  specified  by  any  one  term,  it  may  be  identified  with  that  term, 
which  is  an  affine  function  (an  element  of  A(H)).  Thus,  {constant  sequences}  may  be  identi¬ 
fied  with  the  vector  space  A  ( H ),  and  we  have  the  direct  sum  decomposition: 

PA(r,p)=PA(r,p)0+A(tf). 

One  caution  should  be  given  here.  What  we  have  called  a  polyhedral  surface  may  have 
self-intersections.  It  may  even  flatten  out  into  a  plane,  as  in  the  case  of  a  constant  sequence. 
Under  some  fairly  obvious  conditions,  however,  a  P.A.  sequence  will  give  a  polyhedral  surface 
which  is  a  topological  image  of  the  standard  2-sphere.  For  such  a  "nonsingular"  surface,  we 
may  refer  to  its  interior,  and  we  call  the  closure  of  this  interior  the  corresponding  polyhedron. 
In  our  applications,  we  shall  be  dealing  with  convex  surfaces,  and  these  are  always  non¬ 
singular. 

3.3  The  map  A:  PA(,V,p)-*S  (T,p) 

We  now  define  a  map  A  from  piecewise  affine  sequences  on  T(p)  to  equilibrium 
stresses  on  Hp).  The  map  A:  PA  ,  p)-*S  {T ,  p)  will  be  seen  to  be  a  linear  transforma¬ 
tion.  Recall  that  T  and  p1(  •  •  •  ,p„  are  fixed  Thus,  it  suffices  for  each  etJ  to  express 
u  in  terms  of  the  ar  and  the  vertices  px,  •  •  •  ,p„.  Choose  an  arbitrary  reference  point 
p,  e  H  =  \z  =  1}  that  is  not  collinear  with  any  edge  of  T{p).  Then 

=  e(r,s)(a,(p.)-ar(p.))/[p/,pj,,p»]  .  (11) 

where  [pt.P,,P'J  >8  the  usual  triple  product  in  R3  I f  p,,p;,  and  p.  are  regarded  as  column 
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vectors,  then  [p,  ,p;,p.]  =  det  ( pt,pj,p .)■ 

The  factor  e(r,s)  equals  +  1  or  —  1  according  to  certain  orientation  conventions  Assume 
R3  is  oriented  by  the  "right-hand  rule,"  so  that  when  viewed  *Yom  above,  a  positive  rotation 
in  H  is  counterclockwise  This  counterclockwise  orientation  in  H  imparts  an  orientation  to 
each  interior  face  fr  of  T(p),  hence  an  orientation  to  each  edge  in  the  boundary  of  fr.  We  give 
f0  a  clockwise  orientation. 

Now  consider  the  given  edge  etJ,  i  <j.  It  is  incident  to  precisely  two  faces,  say  fr  and  fs. 
This  defines  the  r  and  s  in  (11).  We  set  e(r,s)  equal  to  + 1  if,  in  the  orientation  imparted  by 
fr,  etJ  points  from  p,  to  p}.  Otherwise,  we  set  e(r,s)  equal  to  —1. 

Note  that  the  value  of  (11)  is  not  affected  by  an  exchange  of  r  and  s,  since  this  changes 
the  sign  of  both  factors  in  the  numerator.  Thus,  the  right-hand  side  of  (11)  depends  only  on  i 
and  j,  so  that  is  well-defined. 

This  concludes  our  description  of  (11). 

Now  (11)  defines  utJ  when  p,  is  adjacent  to  pJt  and  i  < j.  By  using  equations  (1)  (a),(b) 
in  §1,  we  extend  this  to  a  definition  of  utJ  for  1  S  i,  j  tS  n,  i  X  j,  obtaining  a  stress  on 

r 

3.3.1  Theorem:  jw,;|  is  an  equilibrium  stress  on  T(p). 

We  prove  this  in  Appendix  A. 

Formula  (11)  shows  clearly  that  =  A  {ar}  depends  linearly  on  {ar},  i.e.,  that  A  is  a 
linear  transformation  PA  (T.p)  — ►  S  (T,p).  Its  nuilspace  is  easy  to  describe. 

3.3.2  Theorem  The  nuilspace  of  A  consists  of  the  constant  sequences  in  PA(r.p),  i.e.,  using 
our  comments  in  3  2,  nuilspace  A  =  A(H). 

Proof:  Formula  (11)  implies  that  {ar}  is  the  constant  sequence  if  and  onlv  if  the  corr»8Dond- 
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ing  are  all  zero.  f~[ 

Now  using  the  direct  sum  decomposition  in  3  2, 

PA(r,p)  =  PA(r.p)o+A(/f), 

we  get  the  following 

3.3.3  Corollary:  A  I PA(T. p)0 -»  S(T,p)  is  1  -  1.  Q 
Finally,  we  state  the  key  result  of  this  section: 


3.3.4  Theorem:  A:  PA(I\p)  -  S(I\p)  (or  A  |PA(I\ p)0  :  PA(r,p)0  -  SlT.pJi  is  onto. 

The  proof  of  this  theorem  uses  some  elementary  homology  theory;  it  is  presented  in 
Appendix  A. 

Combining  Corollary  3.3  3  and  Theorem  3.3.4,  we  have: 

3.3.5  Corollary:  A  |PA(r,p)0  is  an  isomorphism  PA(I\p)0  -*  SlT.p). 

Corollary  3.3.5  gives  the  1  —  1  correspondence  between  polyhedral  surfaces  and  stressed 
graphs  announced  in  our  introductory  paragraph.  Thus,  schematically,  we  have 


Surfaces  S  over  T(p)  1 

JP.A.  Sequences 

^0 

|. Equilibrium  stresses 

sitting  on  f0 

{{aJs.t.  a0=0 

|  {ul j}  on  f(p ) 

PA(r,p)0  S(T.p) 


where  Ao  =  A  I  PA(r,p)0. 

3.4  Varying  F(p) 

For  our  purposes  we  will  need  to  vary  the  location  of  pi ,  P2.  '  '  ■  ,p„  (the  "geometry"  of 


Dp)),  while  fixing  T  (the  "topology"  of  Hp)).  We  allow  p^  p2,  ■  •  •  ,p„  to  vary  in  H,  sub¬ 
ject  to  the  following  restrictions:. 
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(a)  P1P2.  '  '  •  •  Pn  are  distinct; 

(b)  T(p)  is  a  planar  embedding  of  I\ 

We  modify  the  spaces  in  diagram  (12)  to  make  explicit  the  dependence  on  T  and  p: 


1 13) 


Pairs  (S,  T  (p)) 
where  S  is  a 

Pairs{{ar},  T  (p)) 
s.t.{ar}  is  a  P  A 

Ao 

Pairsifu.j},  T  (p )) 
where  {uLJ}  is  an 

surface  over  T(p) 

sequence  relative  to 

equilibrium  stress 

sitting  on  f q 

T  (p),  with  a0  =  0 

on  f  (p) 

Each  of  these  sets  has  a  natural  topology.  For  example,  suppose  T  has  e  edges.  Then,  each 
pair  ({u^},  F(p ) )  corresponds  to  an  (e  +  3n)-tuple  (  ...  ,  ulJf  ...  ,  px,  ...  ,  p„),  so  that  SfT.p) 
may  be  considered  a  subset  of  R‘  +  3n.  Similarly  for  the  other  sets.  The  left-hand  set  and 
PA(r.p)0  are  essentially  two  versions  of  the  same  thing,  as  described  in  §3.2.  The  left-hand 
arrow  simply  denotes  the  transition  from  one  version  to  another  and  is  thus  clearly  a 
homeomorphism  The  right-hand  arrow  Aq  is  also  a  homeomorphism,  but  this  requires  some 
explanation 

We  would  like  to  define  A0  by  the  formula 


A0(ar,r(p))  =  (wi;,  T(p)),  (14) 

where  {ar}  and  {wIy}  are  related  by  equation  (11)  of  §3.  However,  that  equation  requires  the 
choice  of  a  fixed  reference  point  p»,  and  this  choice  restricts  us  to  p’s  such  that  no  edge  of 
T(p)  is  collinear  with  p«.  Assume  this  restriction  for  the  moment  It  is  then  easy  argue, 
using  3.3.5,  that  A0  gives  a  1-1  correspondence.  Formula  (11)  shows  that  this  correspondence 
and  its  inverse  are  continuous.  That  is,  subject  to  the  above  restriction  on  p,  Ao  is  a 
homeomorphism. 

We  now  argue  that  the  restriction  on  p  may  be  eliminated.  First  note  that  the  restric¬ 
tion  on  p  is  an  "open  condition”;  that  is,  a  specific  choice  of  reference  point  p.  restricts  us  to 
homeomorphic  open  sets  in  PA(r,p)o  and  SlT.p),  say  f/(p»)  and  V(p»),  respectively.  Next 
note  that  by  varying  p«  in  H,  the  resulting  families  {U(p»))  and  { V (p» )}  cover  PA(r,p)o  and 
S(f,p),  respectively  Finally,  note  that  the  value  of  the  right-hand  side  of  (11),  is 
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independent  of  the  choice  of  p«  (see  Appendix  A,  Corollary  A. 5).  This  means  that  the 
homeomorphisms  &q:U (p-) —*  V(p-)  fit  together  compatibly  to  give  a  map 
A(jiPA(r,  p)o  —> >  S(T,p),  which  is,  thus,  a  local  homeomorphism.  But,  using  Corollary  3.3  5 
and  equation  <  14 1  above,  it  is  easy  to  see  that  Aq  is  a  1-1  correspondence.  Thus, 

A0:PA(r,p)0  -S  (r,p) 

is  a  homeomorphism  with  no  restriction  needed  on  p  other  than  that  T(p)  be  a  planar  embed¬ 
ding  of  r 

Sometimes  a  voluntary  restriction  on  p  is  desirable.  For  example,  we  may  wish  to  hold 
some  subset  of  the  vertices  fixed,  while  allowing  the  rest  to  vary.  One  special  case  of  this  is 
important  enough  to  deserve  some  special  mention  This  is  the  case  when  the  fixed  vertices 
comprise  the  peripheral  &-gon  IT  of  T.  Let  PA(r,p)n  and  S(T,p)n  denote  the  closed  sub¬ 
sets  of  PA(r,p)0  and  S(r,p),  respectively,  in  which  the  peripheral  vertices  are  restricted  to 
certain  fixed  values  Then,  A0  in  (14)  restricts  to  a  homeomorphism  (still  denoted  A0) 

PA(r,p)n  tso\p)n. 

For  simplicity,  we  usually  omit  the  subscript  IT  from  the  notation. 

3.5  Remarks  on  stressed  graphs 

3.5.1  It  is  important  for  the  reader  to  be  aware  that  not  every  graph  T(p)  admits  an  equili¬ 
brium  stress  For  example,  the  (solid-line)  graph  in  Figure  3(a)  does  admit  such  a  stress;  the 
(solid-line)  graph  in  Figure  3(b)  does  not. 
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Figure  3.  Examples  of  embedded  graphs  with  and  without  equilibrium  stresses. 

Thus,  whereas  it  is  relatively  easy  to  deform  Tip)  tsay,  from  configuration  * n »  to  tbi 
above),  it  may  not  be  possible  to  deform  an  associated  equilibrium  stress  compatibly  It  is  pre¬ 
cisely  this  phenomenon  that  we  encounter  when  we  attempt  to  deform  polyhedra.  only  now  it 
is  expressed  in  terms  of  stressed  graphs.  In  Section  4  we  show  that  certain  deformations  of 
convex  stressed  graphs  are  possible,  and  this  will  be  sufficient  for  our  applications 

3.5.2  Consider  an  arbitrary  compact,  convex  polyhedron  K  with  face  f q.  By  applying  a  suit¬ 
able  affine  transformation  A,  we  may  arrange  things  so  that  A(/o)  lies  in  the  hyperplane  H 
=  {{x,y,z)|^  =  1}.  We  then  choose  a  point  p  in  R3  outside  the  polyhedron  A  (AT)  but  near 
one  of  the  interior  points  of  A(f0).  By  performing  a  central  projection  from  p,  the  entire 
boundary  of  AIK)  can  be  mapped  into  A(/o),  provided  p  is  close  enough  to  A(f 0)  (see  Figure 
4).  This  is  known  as  a  Schlegel  projection,  and  the  projected  1-skeleton  of  A  (AT)  is  known  aa  a 
Schlegel  diagram  'This  proves  that  such  L-skeieta  are  planar  '  Now,  by  applying  a  projective 
transformation  T  that  fixes  H  and  sends  p  to  the  point  at  infinity  corresponding  to  the  direc¬ 
tion  of  the  r-axis,  we  may  arrange  things  so  that  TA(K)  sits  over  its  face  TA(f o)  in  H 
Moreover  the  Schlegel  projection  is  then  replaced  by  the  usual  orthogonal  projection  onto  H 
Thus,  by  applying  a  suitable  projective  map,  in  this  case  TA,  we  can  transform  any  compact, 
convex  polyhedron  K  and  face  fo  into  the  context  analyzed  in  this  section.  We  shall  call  such 
a  transformation  (or  change  of  coordinates)  a  normalization  of  the  situation  In  our  applica¬ 
tions,  we  usually  deal  only  with  the  normalized  objects  and  leave  it  to  the  reader  to  perform 
the  normalization  or  its  inverse. 
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Figure  4  Schlegel  projection 

3.6  Convex  polyhedra  and  convex  stressed  graphs 

The  following  result  is  well  known  to  experts,  but  a  proof  does  not  seem  to  be  readily 
accessible  in  the  literature.  We  present  a  proof  at  the  end  of  Appendix  A  i  Lemma  A. 9>. 

3.6.1  Theorem:  Suppose  (  T  (p))  is  a  stressed  graph  in  H  and  2  is  the  corresponding 

polyhedral  surface  sitting  over  Tip)  Then  2  is  the  boundary  of  a  convex  polyhedron  if  and 
only  if  Tip)  and  {wi; }  are  convex 

§  4.  Intersecting  a  convex  polyhedron  with  a  half  space 

In  this  section  we  present  an  algorithm  for  intersecting  a  convex  polyhedron  with  a  half 
space  The  algorithm  has  the  property  that  all  conversions  of  numerical  data  to  logical  data 
are  independent  and  consistent.  Moreover,  assuming  a  suitable  degree  of  precision  in  the 
input  data  and  in  the  arithmetic  that  is  used,  it  satisfies  the  definition  of  correctness  of  §1.3 
The  algorithm  is  described  in  §4  2  and  a  proof  of  correctness  is  given  in  §4  4 

§4.1  Representations  and  Accuracy 

A  representation  R  of  a  convex  polyhedron  K  is  a  data  structure  containing  incidence 
information  about  K  as  well  as  face  equations  approximating  those  of  K  We  may  call  K  a 
model  of  R  The  incidence  information,  given  as  logical  data,  lists  the  incidences  among  the 
vertices,  edges,  and  faces  of  K,  what  we  call  the  face -topology  of  K.  The  face  equations,  given 
as  algebraic  and  numerical  data,  correspond  bijectively  to  the  listed  faces,  with  no  a  priori 
assumptions  about  their  degree  of  accuracy  Thus,  a  given  R  may  have  many  models  with 
nothing  in  common  but  their  face-topology 

Topological  criteria  due  to  Tutte  [T]  can  be  used  to  determine  whether  a  candidate  for  a 
representation  actually  is  a  representation  of  some  convex  polyhedron.  Specifically  Tutte 
describes  when  a  given  graph  admits  a  convex  embedding  into  the  plane  with  a  given  sub-k- 
gon  as  periphery  and  possessing  a  restricted,  convex  equilibrium  stress.  Choose  the  periphery 
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to  be  triangular,  if  possible.  Then  the  stress  may  be  modified  to  become  a  convex  equilibrium 
stress  icf.  §2.3 1,  and  the  Maxwell  correspondence  of  §3  implies  that  the  graph  is  the  1-skeleton 
of  a  convex  polyhedron.  If  the  periphery  cannot  be  chosen  to  be  triangular,  the  graph  must 
have  a  degree-three  vertex.  Introduce  a  triangle  in  place  of  this  vertex,  apply  the  foregoing 
argument  to  the  modified  graph,  and  then  restore  the  original  vertex. 

Of  course,  without  some  further  assumptions  on  the  relationship  between  the  incidence 
information  in  R  and  the  face  equations  of  R,  the  notion  of  representation  is  not  very  useful 
For  example,  if  Kq  denotes  the  polyhedron  actually  determined  by  the  face  equations,  Kq  may 
be  the  empty  set,  or  it  may  be  some  unbounded  cone  in  R3,  both  totally  unrelated  to  K. 
Minimally,  we  are  entitled  to  require  that  K0  satisfy  the  following: 

(1)  K o  is  a  (bounded)  convex  polyhedron  with  exactly  one  face  for  each  face  equation  in 
R.  Two  such  faces  meet  only  if  they  are  listed  as  incident  in  R. 

(2)  Each  edge  listed  in  R  actually  appears  in  Kq  as  the  intersection  of  the  appropriate 
faces  Moreover,  the  cyclic  order  of  such  edges  around  a  given  face  occurs  in  Kq  as  indicated 
in  R. 

It  is  easy  to  see  that  if  the  face  equations  of  R  are  sufficiently  close  to  those  of  K,  then 
K o  does  satisfy  (1)  and  (2). 

Note  that  generally  K o  will  have  edges  other  than  those  described  in  >2).  In  fact,  each 
vertex  v  of  R  ior,  equivalently,  of  K)  determines  a  subgraph  T (u)  of  the  1 -skeleton  of  Kq 
Each  edge  of  T (u)  is  the  intersection  of  two  faces  whose  counterparts  in  R  are  not  listed  as 
adjacent  to  each  other  but  are  listed  as  incident  to  v.  It  is  not  hard  to  deduce  from  (1)  and  <2) 
that  the  graphs  T{v)  are  nonempty,  connected,  and  contain  no  cycles.  Thus,  each  T  (v)  is 
either  a  fiingle  vertex  or  a  tree  We  call  it  a  vertex -tree.  By  (1),  every  edge  o(  Kq  belongs  to  a 
vertex-tree  or  arises  as  in  <2). 

The  features  described  in  ill  and  (2)  are  stable  under  small  perturbation  of  face  equa¬ 
tions.  However,  in  general  the  vertex  information  of  K,  carried  by  R,  is  not  preserved  in  Kq 
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because  vertices  of  degree  greater  than  three  are  unstable  with  respect  to  such  perturbation 
It  is  precisely  this  instability  phenomenon  that  has  led  to  most  of  the  considerations  of  this 
paper 

Henceforth  we  consider  only  representations  satisfying  the  two  assumptions  above.  If  we 
wish  to  emphasize  that  they  are  in  force,  we  may  call  the  representation  in  question  a  good 
representation. 

The  algorithm  that  we  present  in  §4  2  will  require  some  additional  considerations  which 
are  more  quantitative  in  nature.  In  general,  the  vertex  trees  T[v)  of  K0  can  be  determined 
exactly  from  R  only  if  infinite-precision  arithmetic  is  used.  In  the  event  that  finite-precision 
arithmetic  is  used,  the  location  of  each  vertex  of  T(u)  may  be  determined  only  up  to  some 
small  radius.  Let  N (v)  denote  the  set  of  all  points  in  R3  whose  distance  from  T (v)  is  less 
than  or  equal  to  this  radius.  To  include  the  infinite-precision  case,  we  allow  the  radiuB  to  be 
zero. 

Now  choose  a  positive  real  number  a.  Then  we  call  R  an  a-representation  of  K  if  R  is  a 
good  representation  and  if  the  following  three  conditions  are  satisfied: 

(3)  For  each  u  in  K,  N (u)  is  contained  in  the  open  ball  B(v,a )  of  radius  a  around  v. 
We  shall  call  B  (v,  a)  a  vertex  ball  belonging  to  v. 

1 4)  No  two  such  balls  meet,  nor  is  any  vertex  ball  with  a  of  the  convex  hull  of  any  two 
others 

(5)  If  a  plane  meets  three  or  more  vertex  balls  belonging  to  vertices  of  a  face,  then  the 
only  vertex  balls  it  meets  belong  to  vertices  of  that  face 

Condition  (3)  gives  a  measure  of  the  accuracy  of  the  representation  with  respect  to  the 
model  K.  It  asserts  that  the  total  imprecision  coming  from  both  the  face-equation  data  of  R 
and  the  finite-precision  computations  that  may  be  used  is  less  than  a.  Condition  >4)  asserts 
that  the  edges  and  faces  of  K  are  large  relative  to  a  (or  inversely,  that  a  is  small  relative  to 
these  features).  It  also  implies  that  the  angles  between  consecutive  edges  in  a  face  are 
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bounded  away  from  it  or  0  and  that  the  dihedral  angle  between  a  face  and  a  plane  passing 
near  three  of  its  vertices  is  bounded  away  from  it' 2.  Condition  <5)  relates  a  to  the  dihedral 
angles  of  K  It  will  be  violated  if  and  only  if  these  angles  are  too  close  to  0  or  it  relative  to  a 

Given  K,  conditions  (4)  and  (5i  may  be  formulated  as  upper  bounds  on  a.  Thus,  for  each 
convex  polyhedron  K,  there  exists  a  positive  number  (i  such  that  K  has  a-representations  for 
all  a  <  fi.  If  if  is  a  representation  of  K  and  the  total  imprecision  is  known  a  prion,  then  the 
upper  bounds  of  i4)  and  (5)  may  be  estimated  from  R.  When  both  estimates  are  greater  than 
the  total  imprecision,  then,  for  any  choice  of  a  in  between,  R  is  an  a-representation.  If  one  of 
the  estimates  is  less  than  the  total  imprecision,  then  we  must  reject  the  representation  as 
being  insufficiently  precise. 

$4.2  The  algorithm 

Let  a  be  a  small  positive  number.  The  algorithm  takes  as  input  an  a-representation  R 
of  a  convex  polyhedron  K  and  a  linear  inequality  X  2t  0  defining  a  half  space.  Its  output  is  a 
representation  R*  of  a  convex  polyhedron  K* .  We  may  call  R*  the  intersection -representation. 

The  algorithm  proceeds  as  follows.  Let  P  be  the  plane  X  =0.  For  each  vertex  v  of  K, 
the  algorithm  attempts  to  determine  by  means  of  a  numerical  computation  the  side  of  P  on 
which  v  lies.  Such  a  determination  will  succeed  exactly  when  the  computed  location  of  v  has 
distance  from  P  greater  than  or  equal  to  a  in  this  case,  the  algorithm  reports  the  side  con¬ 
taining  the  computed  location.  If  the  attempted  determination  fails,  that  is,  if  the  computed 
location  is  strictly  within  a  of  P,  then  u  is  reported  to  be  on  P 

If  none  of  the  vertices  of  K  are  reported  on  the  side  X  >  0,  then  the  algorithm  returns 
the  null  representation.  Similarly  if  none  are  on  X  <  0,  then  the  algorithm  returns  R. 

Suppose  that  three  or  more  vertices  of  a  face  F  are  reported  on  P,  and  some  vertices  of  K 
are  reported  on  each  side  of  P.  By  property  (5)  of  an  a-representation,  this  can  occur  for  at 
most  one  face  F  of  K.  In  this  case,  the  algorithm  computes  the  inward  pointing  normals  of  F 
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and  P  (i.e  ,  for  F  this  means  pointing  inward  into  K,  for  P  this  means  pointing  inward  into 
A  ^  0),  and  then  it  computes  the  dot  product  of  these  normals.  It  returns  the  null  representa¬ 
tion  if  the  dot  product  is  negative  and  R  if  it  is  positive.  Note  that,  by  property  <4)  of  an  a- 
representation,  the  absolute  value  of  this  dot  product  is  bounded  below.  Therefore,  with  only  a 
modest  assumption  on  total  precision,  that  is,  on  the  input  a,  the  algorithm  is  able  to  make 
the  decision  in  this  case. 

So  let  us  now  assume  that  each  face  has  at  most  two  vertices  reported  on  P  and  that 
some  vertices  of  K  are  reported  on  each  side  of  P.  To  construct  the  incidence  data  for  the 
intersection-representation  in  this  case,  partition  each  edge  of  K  that  intersects  P.  This  is 
done  solely  as  a  logical  computation  based  on  the  location  of  endpoints  of  the  edge  relative  to 
P.  Discard  segments  on  the  "wrong"  side  of  P,  and  to  the  remaining  set  of  new  edges,  adjoin 
all  the  edges  of  K  both  of  whose  endpoints  are  reported  to  be  on  P  or  on  the  correct  side  of  P. 
Finally,  determine  new  edges  by  joining  appropriate  partition  points  constructed  above,  or  a 
partition  point  with  a  vertex  reported  on  P,  as  the  case  may  be.  ThiB  gives  all  the  edges  and 
vertices  for  the  intersection-representation.  The  faces  are  of  three  kinds:  the  face  determined 
by  P  itself,  portions  of  faces  of  K  which  are  divided  by  P,  and  faces  of  K  lying  on  the  "correct" 
side  of  P.  Again,  these  are  all  determined  by  logical  computation,  as  are  the  required 
incidences.  For  the  face  determined  by  P,  we  have  the  face  equation  A  =  0.  For  the  remain¬ 
ing  faces,  the  face  equations  come  directly  from  R 

This  concludes  our  description  of  the  algorithm  and  of  the  output  representation  R*.  We 
need  to  show,  however,  that  R *  actually  is  a  representation. 

4.2.1  Theorem;  R*  is  a  good  representation. 

Proof:  The  null  representation  and  R  itself  are  both  good  representations,  so  we  may  assume 
that  the  algorithm  has  reported  vertices  on  both  sides  of  P,  at  most  two  of  which  are  on  a 
given  face 
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Recall  that  Kq  is  the  bounded,  convex  polyhedron  given  by  the  face  equations  listed  in 
R  Let  K*  o  be  the  analogous  convex  polyhedron  for  R* .  To  show  that  R *  is  a  representation, 
we  refer  the  reader  to  the  second  paragraph  of  §4.1  and  the  result  of  Tutte  mentioned  there. 
Tutte’s  result  asserts  that  it  is  sufficient  to  show  that  the  edge-graph  G*  given  in  R*  is  3- 
connected  and  planar.  Let  G*  q  denote  the  edge-graph  of  K*0,  by  Tutte  [T],  together  with  the 
Maxwell  correspondence  of  Section  3,  G*  q  has  these  properties.  Moreover,  G*  may  be 
obtained  from  G*  0  by  collapsing  each  vertex  tree  to  a  point  .  It  is  not  hard  to  show  that  col¬ 
lapsing  an  edge  of  a  vertex  tree  does  not  destroy  planarity  (e.g.,  the  reverse  operation  cannot 
eliminate  Kuratowski  subgraphs).  Thus  G*  is  planar.  To  prove  that  it  is  3-connected,  observe 
that  the  link  of  every  vertex-tree  is  a  simple  closed  curve.  (By  the  link  of  a  subset  S  of  G*0, 
here,  we  mean  the  union  of  all  edges  e  of  G*  o  such  that  e  is  disjoint  from  S  but  e  belongs  to  a 
face  meeting  5.)  Collapse  a  single  tree  to  a  point  p\  this  does  not  affect  the  link.  If  p  and 
another  point  are  selected  in  the  resulting  graph,  the  3-connectivity  of  G*  0  implies  that  each 
component  of  the  complement  of  these  points  meets  the  link  of  p.  Since  one  point  cannot 
disconnect  a  simple  closed  curve,  this  proves  that  the  complement  is  connected.  This  same 
conclusion  follows  even  more  easily  for  the  complement  of  two  points  each  different  from  p. 
Thus,  the  collapse  of  a  single  tree  in  G*  o  results  in  a  3-connected  planar  graph.  Now  con¬ 
tinue  this  process,  collapsing  all  the  trees  in  G%. 

It  remains  to  show  that  R  satisfies  conditions  (1)  and  (2)  that  define  a  good  representa¬ 
tion  We  first  observe  that  K*  o  is  bounded.  To  see  this,  consider  the  subset  M  of  Kq  on 
which  the  linear  function  \  assumes  its  maximal  value.  M  is  either  a  single  point,  a  line  or  a 
single  face  of  Kq.  In  any  case,  M,  P,  and  the  planes  of  faces  of  Kq  incident  to  M  enclose  a 
bounded  region  that  contains  K*q.  Clearly  every  face  equation  listed  in  R*  gives  a  face  of 
K*  o,  because  such  an  equation  either  corresponds  to  a  face  or  portion  of  a  face  of  Kq  (since  R 
is  good),  or  it  corresponds  to  the  face  determined  by  P.  By  construction,  K*  0  has  no  addi¬ 
tional  faces.  If  two  faces  of  K*  q  meet  in  an  edge,  then  by  the  definition  of  R*  and  the  fact 
that  R  is  good,  these  faces  are  listed  as  incident  in  R*.  Thus  R*  and  K* q  satisfy  property  (1) 
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of  good  representations.  For  property  (2),  consider  the  edges  listed  in  R*  that  correspond  to 
edges  of  K,  or  portions  of  edges  of  K,  but  do  not  have  both  endpoints  reported  on  P  Such 
edges  occur  in  K0  as  intersections  of  appropriate  faces  and  hence  occur  in  K*  0.  The  remain¬ 
ing  edges  of  R*  are  portions  of  intersections  of  face  planes  of  Kq  with  P  That  these  edges  all 
occur  in  the  appropriate  order  in  the  boundaries  of  faces  follows  from  the  corresponding  fact 
for  K o  and  the  construction  of  R* .  Thus,  R *  and  K*  o  satisfy  (2).  Q 

By  Theorem  4.2.1,  the  algorithm  produces  a  valid,  consistent  geometric  object.  It 
remains  to  demonstrate  that  our  algorithm  is  correct  for  the  operation  of  intersection,  as  per 
the  definition  of  §1.3.  In  §4.4,  we  show  that  this  is  the  case,  provided  that  a  is  sufficiently 
small.  It  is  also  desirable  to  improve  Theorem  4.2.1  so  that  we  know  that  R*  is  not  merely  a 
good  representation  but  also  that  it  is  a  /3-representation,  where  /3  is  suitably  related  to  a. 
We  also  do  this  in  §4.4  again  under  the  assumption  that  a  is  sufficiently  small.  In  §4.3  we 
prove  some  of  the  technical  results  necessary  for  these  demonstrations. 

4.3  The  main  technical  results 

In  this  section,  we  present  a  technical  result.  Theorem  4.3.3,  which  we  need  for  our 
demonstration  that  the  algorithm  of  §4.2  is  correct.  Actually,  the  demonstration  requires  a 
similar  but  more  general  result,  which,  because  its  proof  is  substantially  more  complicated,  we 
present  in  Appendix  C. 

We  continue  to  use  the  terminology  and  notation  of  §4.2.  Thus,  in  particular,  K  denotes 
the  convex  polyhedron  of  the  algorithm,  and  X  denotes  the  given  linear  form.  We  cannot 
move  vertices  of  K  arbitrarily,  even  by  small  amounts,  and  expect  to  find  another  polyhedron 
with  those  vertex  locations  and  the  same  face-topology  as  K.  However,  Theorem  4.3.3  tells  us 
that  we  can  make  a  small  arbitrary  vertex  change  in  the  "X-direction"  without  disturbing  the 
face-topology  of  K,  provided  we  allow  an  accompanying  change  in  other  directions. 

The  proof  of  Theorem  4.3.3  deals  entirely  with  a  stressed  equilibrium  graph  for  K.  Thus, 
we  suppose  that  we  have  normalized  the  situation  as  discussed  in  §3.5.2.  In  particular,  K  sits 
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over  a  convex  3-connected  graph  T(p)  in  the  plane  H  =  {(x,y,x)|z  =  1},  and  P  is  the  plane 
X  =  0.  We  suppose,  as  usual,  that  the  vertices  p,  of  the  graph  have  been  numbered  so  that 
P1P2-  '  '  -P*  ar«  peripheral  and  the  rest  interior.  Call  an  edge  of  Hp)  vertical  if  its  end¬ 

points  have  the  same  x-coordinate  The  proof  of  Theorem  4  3.3  makes  use  of  the  fact  that 
every  such  graph  contains  a  spanning  tree  with  no  vertical  edges.  We  prove  a  slightly 
stronger  result  in  Lemma  4.3.2,  but  first  we  show  how  such  a  spanning  tree  is  used. 

4.3.1  Lemma:  Let  T(p)  be  as  above  and  suppose  that  it  has  a  spanning  tree  T  with  no  verti¬ 
cal  edges.  For  each  edge  etJ  nut  in  T,  let  be  an  arbitrary  stress  constant.  Then  there 
exists  a  unique  assignment  of  stress  constants  to  the  edges  of  T,  hence  a  stress  on  T, 
such  that  the  x-equilibrium  equation 

2  «,,u,  ~xj]  =  0 

j =1 

is  satisfied  for  each  1  =  1,  2,  ....  n.  Here,  x,  denotes  the  x-coordinate  of  p,. 

Proof:  Suppose  that  pi  is  a  leaf  of  T  and  t?i2  is  an  edge  of  T.  All  the  stress  constants  «ly, 
j  >  2,  are  already  determined,  and  Xi  —  X2  *  0.  Thus,  the  above  equation,  for  1  =  1,  can 
be  used  to  solve  for  The  remaining  stress  constants  are  found  similarly  starting  at  the 

leaves  of  T  and  working  towards  the  root.  Then,  x-equilibrium  is  guaranteed  by  construction 
at  every  vertex  except  possibly  at  the  root  of  T.  But,  in  fact,  x-equilibrium  holds  there  as 
well,  because  the  sum  of  all  the  expressions 

n 

2  -*,) 

j  =  1 

i  =  l,2 . n,  is  formally  equal  to  zero.  □ 

Next  we  show  that  not  only  does  T  (p)  contain  a  spanning  tree  with  no  vertical  edges 
but  it  also  contains  one  with  no  edges  that  are  near-vertical.  We  say  that  an  edge  is 
near  —  vertical  if  its  endpoints  have  x-coordinates  that  differ  by  less  than  some  prescribed, 
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small,  positive  quantity  v.  The  precise  choice  of  v  is  not  important.  We  require  only  that  no 
face  of  T  (pJ  have  consecutive  edges  that  are  near-vertical.  Note  that,  by  the  convexity  of  K 
and  the  fact  that  it  sits  over’  T  (p),  there  exist  sufficiently  small  v  such  that  this  property  is 
satisfied. 

4.3.2  Lemma:  Let  T  (p)  be  as  above,  and  assume  that  at  most  one  of  its  peripheral  edges  is 
near-vertical.  Then  T  (p)  has  a  spanning  tree  with  no  near-vertical  edges 

Proof:  We  need  only  show  that  the  deletion  of  ail  near- vertical  edges  leaves  T(p)  connected. 
Since  at  most  one  peripheral  edge  is  near-vertical,  the  modified  periphery  remains  connected. 
Since  no  vertex  is  incident  to  a  pair  of  consecutive  near- vertical  edges  and  T  (p)  is  convex  and 
3-connected,  each  interior  vertex  is  incident  to  some  non-near-vertical  edge  going  to  the  left. 
Thus,  each  interior  vertex  may  be  connected  to  the  modified  periphery  by  a  leftward  moving 
path  consisting  of  non-near-vertical  edges.  Q 

Note  that  when  the  periphery  of  T  (p)  is  triangular,  then  it  automatically  satisfies  the 
hypothesis  of  Lemma  4.3.2. 

4.3.3  Theorem:  Let  K  be  as  above,  and  let  Xo  denote  the  n-tuple  of  x-coordinates  of  the  ver¬ 
tices  of  K.  There  exists  an  5  >  0  such  that  if  X  j  is  any  n-tuple  satisfying 

|X0  -  Xi  |  <  8, 

then  there  ia  a  convex  polyhedron  K  j  with  the  same  face-topology  as  K  and  having  X  i  as  the 
n-tuple  of  x-coordinates  of  its  vertices.  Furthermore,  K  i  may  be  chosen  so  that  its  vertices 
vary  smoothly  as  a  function  of  X  i . 

Proof:  Let  ( {u;j  /,  T  (pJ  )  be  the  convex  stressed  graph  in  H  corresponding  to  K  For  simpli¬ 
city,  we  deal  here  only  with  the  case  in  which  the  periphery  of  T  (p)  is  a  triangle.  A  more 
general  situation  is  considered  in  Appendix  C. 
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Set  A*o  —  U  i  .  x2  •  -.-tr,  )  and  At  —  (x  i  ,  x'2  .  ...,x'„  ).  Since  the  periphery  of  T  (p) 
is  non-degenerate,  the  same  will  be  true  when  xj  ,x2  x3  are  replaced  by  x \  .x  2  ,x'3, 
respectively,  provided  |Ao  —  A’i  |  is  suitably  small.  Our  first  restriction  on  Aj,  then,  is 
that  it  satisfy  |  A’ q  —  A'  i  |  <  6  ^ ,  for  such  a  suitably  small  . 

Next,  choose  a  spanning  tree  T  with  no  near-vertical  edges,  and  let  fi2  =  min 

/lxi  ~  xj  l:  eij  ‘s  an  edge  of  T  }.  Our  second  restriction  on  A'!  is  that 
|A0  ~  K\  |  <  52  2.  This  insures  that  no  edge  of  T  becomes  vertical  when  the  x- 
coordinates  of  its  endpoints  are  replaced  by  the  corresponding  coordinates  of  A  i 

We  then  define  a  new  stress  {u ' LJ  }  for  T  as  follows.  For  every  edge  e,7  not  in  T,  set 
'  ij  —  u)tJ.  For  the  remaining  edges,  proceed  as  in  Lemma  4.3.1,  using  T  and  the  coordinates 

of  Aj.  Clearly  the  resulting  stress  is  a  smooth  function  of  At.  Since  /cj',7  /  =  (utJ  }  when 

A" !  =  A0,  and  (u,j  /  is  a  convex  stress  (i.e.,  utJ  is  negative  for  etJ  peripheral  and  positive 
otherwise),  there  exists  a  positive  63  such  that  /  is  a  convex  stress  for  all  A3  satisfying 
|A0  —  Ai  |  <  83.  Our  final  restriction  on  A !,  then,  is  that  it  satisfy  this  last  inequality. 

In  summary,  then,  the  positive  quantity  5  of  the  theorem  may  be  defined  as 

min^Si  ,  82  }■ 

Now  set  p ' r  —  (x'r  ,yr  ,  1 )  e  H,  r  =  1,2,3,  where  yr  is  the  y-coordinate  of  pr.  By  the 
choice  of  8,  these  are  the  vertices  of  a  non-degenerate  triangle  in  H.  Since  {u }  is  a  convex 
stress  on  T,  we  may  apply  Theorem  2.1,3  to  obtain  unique  points  p'r  ,r  >  3,  satisfying:  (1) 
{u  ’1;  /  is  a  restricted  convex  equilibrium  stress  for  T  (p  );  (2)  the  p  r  depend  smoothly  on  the 
u  tJ  .  1  >  3  or  j  >3  (hence,  smoothly  on  A^.  We  claim  that  Aj  is  precisely  the  n-tuple  of 
x-coordinates  of  p  .  For,  by  construction,  the  coordinates  of  A3  satisfy  x-equilibrium  with 
respect  to  {u ' /  In  matrix  form,  this  is  just  the  first  of  the  two  equations  (8)  of  §2,  or, 
equivalently,  the  first  of  the  two  equations  (9>.  But,  by  Corollary  2.2.2,  this  same  equation 
gives  the  x-coordinates  of  p  4  ,p  5  ,  .  .  .  ,p in  terms  of  the  x-coordinates  of  p' \  ,p  2  ,p  3  , 
wt  ich  coincide  with  x'i  ,x'2  .  *  3  .  respectively,  by  definition.  Thus  the  claim  is  verified. 


Now  modify  the  stress  constants  u  kJ  ,i,j  <  4,  as  described  in  Theorem  2.3  1.  so  as  to  obtain 
a  convex  equilibrium  stress,  still  denoted  (u  tJ  }.  on  T  ip').  By  (2>  above,  and  by  assertion  1 2 ) 
of  Theorem  2.3.1,  the  stressed  graph  ({u‘ LJ  },  T(p '))  depends  smoothly  on  Xj. 

The  results  of  S3  then  produce  a  convex  polyhedron  K  i  corresponding  to  ({u  tJ  /,  F(p)). 
The  last  assertion  of  the  previous  paragraph,  together  with  the  results  of  §3,  shows  that  the 
vertices  of  K\  depend  smoothly  on  X|.  Moreover,  they  have  the  correct  x-coordinates,  by  con¬ 
struction.  0 

4.3.4  Remark:  Since  the  vertices  of  depend  smoothly  on  Xi,  it  is  clear  that  we  may 
choose  K i  to  be  close  to  K  by  choosing  Xi  suitably  close  to  Xo.  More  precisely,  for  every 
t  >  0,  there  is  a  5  >0  such  that  if  |X0  —X!  |  <5,  then  vertices  of  K\  are  within  e,  of 
the  corresponding  vertices  of  K.  To  compute  such  a  5  in  terms  of  e,  one  must  factor  the 

smooth  map  Xx - *KX  into  its  three  component  parts, 

Xi - \j  /,  {u‘ } - *p  ,  and  (/u'M  },  Tip  )) - tXj,  and  then  estimate  moduli  of 

continuity  for  these.  For  our  purposes  here,  it  is  sufficient  to  note  that,  given  K  and  e,  these 
steps  can  be  carried  out.  In  addition,  one  can  factor  in  the  effect  of  normalization  and  its 
inverse.  A  similar  but  more  complicated  analysis  applies  to  the  more  general  method  used  in 
Appendix  C  In  §4  4  we  discuss  how  8  may  be  estimated  when  only  R  is  known. 

4.4  Correctness  and  Accuracy 

To  establish  the  correctness  of  our  intersection  algorithm  according  to  the  definition  of 
§1.3,  it  is  sufficient  to  show  that  a  model  for  the  input  a-representation  R  exists  whose  actual 
intersection  with  the  half-space  A^O  is  a  model  for  the  output  representation  R *.  Note  that 
we  are  not  imposing  any  a  priori  standard  of  accuracy  on  R*  Of  course,  we  could  do  so  by 
requiring  that  R*  be  a  fi  —  representation  of  the  actual  intersection,  for  a  prescribed  >  0. 
We  touch  on  this  below.  In  any  case,  we  are  left  with  the  problem  of  producing  a  suitable 
model  for  the  input  representation.  In  those  cases  in  which  the  algorithm  produces  the  null 
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representation  or  R  itself,  this  problem  is  easily  solved.  For  in  those  cases,  either  a  small 
translation  will  perturb  K  so  that  it  no  longer  meets  the  plane  P  or  small  rotation  followed  by 
translation  has  this  effect.  In  either  case,  this  perturbation  of  K  is  then  the  desired  model  of 
R 

For  the  generic  case  when  vertices  are  reported  on  both  sides  of  P  and  no  face  has  more 
than  two  vertices  reported  on  P,  we  again  start  with  the  given  polyhedron  K  and  attempt  to 
perturb  it  slightly  so  as  to  meet  our  requirements  That  is,  we  attempt  to  perturb  K,  without 
changing  its  face-topology,  so  that  precisely  those  vertices  that  are  reported  on  P  are  actually 
moved  on  to  P  and  no  vertex  moves  from  one  side  of  P  to  the  other.  If  we  can  do  this,  then 
the  perturbed  image  of  K  is  again  clearly  the  desired  input  model. 

Let  us  say  that  a  face  is  reported  to  meet  P  if  it  has  vertices  reported  on  both  sides  of  P. 
Suppose  first  that  a  triangular  face  F  is  reported  to  meet  P  We  then  normalize  the  situation 
as  in  §4.3,  projecting  K  onto  F,  moving  F  to  the  plane  H  and  P  to  the  plane  x  =  0.  The  pro¬ 
jected  image  of  P  in  H  is  a  line  L  that  divides  the  graph  T  (p).  Vertices  within  a  of  P  are 
projected  to  vertices  in  T(p)  within  a  of  L.  We  attempt  to  move  these  vertices  onto  L ;  that  is, 
to  change  their  x-coordinates  to  0  By  Theorem  4.3.3,  this  will  be  possible  if  a  is  less  than  a 
certain  quantity  8 

A  similar  argument  may  be  used  if  P  fails  to  meet  a  triangular  face  as  above  but  is  suit¬ 
ably  close  to  an  order-three  vertex.  For  then  we  need  only  slice  off  this  vertex  in  a  direction 
transversal  to  P,  then  apply  the  preceding  argument,  and  finally  restore  the  order-three  ver¬ 
tex. 

We  next  consider  the  remaining  possibilities  It  is  possible,  for  example,  that  there  are 
faces  reported  to  meet  P  none  of  which  is  triangular  and  P  is  not  close  to  an  order-three  ver¬ 
tex.  In  this  case,  the  same  discussion  applies  as  in  the  previous  paragraph,  only  we  use  the 
more  general  result  of  Appendix  C  instead  of  Theorem  4.3.3.  The  same  conclusion  holds. 
Finally,  if  no  face  is  reported  to  meet  P  and  no  order-three  vertex  is  close  P.  but  there  are  ver- 
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tices 

reported  on  both  sides  of  P,  then  K  must  be  divided  onto  two  pieces  by  a  simple  closed  curve 
consisting  of  edges  that  lie  on  P  or  nearly  lie  on  P  (i.e.,  each  endpoint  of  such  an  edge  is 
reported  to  be  on  Pi.  In  this  case,  we  choose  such  an  edge  and  perturb  K  slightly  so  that  the 
edge  is  moved  onto  P.  Then  we  choose  either  of  the  two  faces  f  \  and  fi  containing  the  edge 
and  project  K  onto  the  plane  of  this  face,  say  f\,  using  a  central  projection  from  a  point  close 
to  f i,  not  in  K  but  on  the  plane  P  (as  usual).  Normalize  as  before  so  that  the  plane  of  f  j 
becomes  H,  P  becomes  x  =  0,  and  the  projection  is  orthogonal.  We  get  a  stressed  graph  T  (p) 
as  before,  only  its  periphery  now  encloses  the  image  of  both  f\  and  f%  which  are  exterior  to 
one  another.  The  faces  are  all  convex,  the  interior  stress  constants  are  positive,  and  the  peri¬ 
pheral  ones  are  negative,  just  as  before.  This  case  can  be  handled  just  as  the  previous  one  and 
is  discussed  in  more  detail  in  Appendix  C. 

Thus,  in  all  cases,  we  can  obtain  the  desired  deformation  provided  a  is  small  enough. 
We  now  formulate  this  as  a  theorem. 

4.4.1  Theorem:  Given  any  convex  polyhedron  K  and  linear  form  X,  there  is  a  8  >  0  such 
that  the  intersection  algorithm  of  §4.2  is  correct  for  all  a-representations  R  of  K  satisfying 
a  <  8  Moreover,  8  may  be  chosen  to  be  independent  of  X. 

Proof:  Only  the  second  assertion  still  requires  verification.  Without  loss  of  generality,  we 
restrict  attention  to  forms  X  =  ax  +  by  +  cz  +  d,  such  that  a2+62  +  c2  =  l.  Those 
X  for  which  the  output  of  the  algorithm  is  the  null  representation  or  R  require  no  6:  that  is, 
no  restriction  on  a  is  needed  for  those  branches  of  the  algorithm  The  remaining  X  satisfy  the 
condition  that  vertices  are  reported  on  both  sides  of  X  =  0  but  at  most  two  from  each  face  are 
reported  on  X  =  0  This  is  a  closed,  bounded  condition  on  the  coefficients  of  X,  and  so  this  set 
of  X  comprise  a  compact  set  C. 
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Choose  a  particular  At  in  C,  and  let  fix  be  the  value  of  6  given  by  the  discussion  above, 
i.e.,  by  Theorem  4.3.3  after  normalizing  the  situation.  Let  us  say  that  a  value  of  6  is  suitable 
for  A  if,  whenever  a  <  S,  K  may  be  perturbed  so  that  precisely  the  vertices  within  a  of 
A  =  0  are  moved  onto  that  plane  and  no  verter  changes  sides.  We  claim  that  fit  2  is  a  suit¬ 
able  value  of  5  for  all  A  in  some  neighborhood  of  At  in  C.  To  verify  this  claim,  choose  an 
a  <  fix'2  and  select  those  vertices  of  K  within  a  of  A  =  0.  Let  T(A)  denote  the  projective 
transformation  that  we  use  to  normalize  K.  For  A  close  to  At  ,  T( A)  may  be  chosen  to  be  a 
smooth  function  of  A,  and  hence  T( A)  is  close  to  TlAi).  Because  the  x-coordinates  of  the 
selected  vertices  of  T(\)(K)  are  less  than  a,  it  then  follows  that  the  x-coordinates  of  the 
corresponding  vertices  of  T(\i)(K)  are  less  than  S\.  Therefore,  as  in  Theorem  4 .3.3,  there 
exists  a  perturbation  h  of  T(\\)(K)  sending  the  selected  vertices  onto  x  =  0,  etc.  Now  when 
A  is  close  to  Aj,  there  is  an  affine  transformation  of  R3  that  preserves  orientation  and 
transforms  Aj  to  A:  A  =  Ax  A.  Define  a  mapping  of  g  of  R 3  by 

g  —  A-1  T(Ax  )_1  /iT(Ax),  an<*  consider  any  of  the  selected  vertices  p,  of  K.  Since 
TTAiHp,)  is  taken  onto  x  =  0  hv  h,  and  T( Ax)  takes  Ax  =  0  to  x  =  0,  \g(pt)  =  0,  as 
desired.  Similarly  g  preserves  the  half-spaces  determined  by  A  =  0.  Thus,  g(K)  is  a  pertur¬ 
bation  of  K  with  desired  properties,  which  verifies  the  claim. 

Since  finitely  many  such  neighborhoods  cover  C,  each  with  a  suitable  value  of  fi,  the 
smallest  of  these  values  is  suitable  for  all  A  in  C,  and  therefore,  for  all  A.  lj 

We  now  turn  to  the  question  of  the  accuracy  of  the  representation  produced  by  the  algo¬ 
rithm.  There  are  (at  least)  two  senses  in  which  the  term  accuracy  can  be  used  here  First,  we 
can  ask  how  closely  the  output  representation  R*  represents  some  model.  Secondly,  we  can 
ask  how  closely  such  a  model  approximates  the  true  intersection  of  K  and  A  2  0. 

The  first  question  has  a  relatively  straightforward  mathematical  answer 

4.4.2  Theorem:  There  exists  a  positive  constant  k  depending  only  on  K,  such  that  if  R  is  an 
a -representation  of  K  and  a  is  less  than  the  constant  5  of  Theorem  4.4.1.  then  each  vertex  of 
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the  model  K *  of  ft  *  constructed  as  in  the  correctness  proof  above  is  within  ica  of  its  location 
computed  using  R " 

Proof:  The  proof  consists  of  looking  closely  at  the  construction  of  the  modified  input  model  in 
our  demonstration  of  the  correctness  of  the  algorithm.  Suppose  first  that  the  algorithm 
reports  all  vertices  on  one  side  of  P  or  on  ft.  If  this  side  corresponds  to  X  2  0,  then  K*  equals 
a  small  translate  of  K  into  this  side  of  P.  The  length  of  the  translation  need  be  at  most  2a,  so 
that  R*  =  ft  is  a  2a-representation  of  K*  in  this  case.  If  the  side  reported  corresponds  to 
A  2  0,  then  K*  equals  the  empty  polyhedron. 

The  case  in  which  the  algorithm  reports  at  least  three  vertices  of  a  face  ft  on  ft  and  some 
vertices  on  both  sides  of  ft  is  more  complicated.  In  this  case,  we  rc'ate  K  through  a  small 
angle  around  the  intersection  line  ft  H  ft  until  ft  lies  on  ft,  and  then  we  perform  a  small 
translation,  say  of  length  a,  making  the  polyhedron  disjoint  from  ft  If  the  result  lies  in 
X  2  0  then  K *  equals  this  resulting  polyedron.  Otherwise,  K*  is  empty.  It  remains  to  show 
that  there  is  a  constant  M  such  that  each  vertex  of  K*  is  less  than  M  a  from  the  correspond¬ 
ing  vertex  of  K.  For  then  the  distance  between  such  a  vertex  and  that  computed  via  ft*(  =  ft) 
is  less  than  (M  +  l)a,  as  desired.  But,  for  any  fixed  set  of  three  or  more  vertices  of  ft  within 
a  of  ft  and  line  segment  of  intersection  of  ft  and  ft.  the  distance  moved  by  vertices  under  the 
rotation  is  a  smooth  function  of  the  angle  of  rotation,  and  this  angle  in  turn  has  an  upper 
bound  of  the  form  m  a  Similarly,  the  locations  of  the  endpoints  of  the  line  segment  are  deter¬ 
mined  up  to  some  m  a.  Thus,  there  is  a  constant  M'  such  that  M  a  bounds  the  distance 
moved  by  vertices  under  rotation  An  additional  a  takes  care  of  the  small  translation,  so  that 
M  =  M'  +  1  is  the  desired  constant.  Note  that  we  may  repeat  this  argument  for  any  set  of 
three  or  more  vertices  on  any  face,  thus  finding  a  single  M  independent  of  X  in  this  case. 

Next  we  choose  any  X  in  the  compact  set  C  described  in  the  proof  of  the  previous 
theorem.  Let  T  be  the  projective  transformation  that  maps  K  into  its  normalized  position  over 
the  graph  f  (p)  in  H  For  any  polyhedron  K  with  n  vertices,  let  u(K')  denote  the  n-tuple  of 
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vertices,  and  let  x{K  )  be  the  n-tuple  of  x-coordinates  of  these  vertices  Let  r  be  the  modulus 
of  continuity  of  T  on  some  large  closed,  bounded  region  Now  apply  the  construction  of 
Theorem  4  3  3  or  Appendix  C,  as  the  case  may  be  to  the  polyhedron  T(K)  and  some  n-tuple 
of  x-coordinates  u  =  ( u  \  ,  u  o  ,  .  .  .  ,  un  ).  The  result  is  a  polyhedron  A"  i  satisfying 
x(A'  i  )  =  u  and  depending  smoothly  on  u.  Thus,  for  u  ranging  in  a  compact  ball  around 
x(T(K)).  there  is  a  constant  .V/  such  that  |r(7’(AO)  —  u(K^  )|  <  M\x(T(K))  —  u  |.  In 
our  perturbation  of  A'  to  construct  a  model  of  R* ,  the  vertices  that  are  moved  onto  P  move  no 
more  than  a  units  in  the  direction  orthogonal  to  P  Since  T  is  just  a  rigid  motion  followed  by 
a  central  projection  from  a  point  close  to  K  but  outside  of  K,  the  coordinates  in  x(T(K))  of 
vertices  that  we  want  to  perturb  are  all  bounded  by  a.  For  our  construction  of  the  desired 
model,  we  set  u  =  0.  Thus,  we  get  jt'(T(K))  —  v(K\  )|  <  M  a.  Next,  we  apply  the 
inverse  transformation  P_1  to  T{K)  and  K\  and  conclude  that 
|  v\K)  —  v(T  (Ki  ))  j  <  .Wt'1  a.  Since  R  is  an  a-representation  of  K,  we  see  that  a 
vertex  computed  using  R  is  within  (Mr-1  4-  Da  of  the  corresponding  vertex  in  T~l  (K \  ). 
Finally,  recall  that  to  construct  the  model  for  the  output  R*,  we  may  also  have  to  introduce 
new  vertices  where  edges  of  T  ~  1  {K\  )  cut  the  plane  P.  Such  an  edge  has  endpoints  at  least 
pa  units  from  P,  for  some  p  that  depends  only  on  T,  and  so  the  angle  it  forms  with  P  is 
bounded  below  It  follows  that  the  indeterminacy  for  such  vertices  is  bounded  by  N  a  for  some 
constant  N  Thus,  setting  K*  equal  to  the  actual  intersection  of  A  2  0  with  T~l  (K\  )  and 
K  =  max/.V,  \f  r  ~ 1  +  1/,  we  see  that  the  vertices  of  K*  are  within  it  a  of  those  computed 
using  R  *. 

We  can  now  combine  the  foregoing  argument  with  a  compactness  argument  similar  to 
that  ot  Theorem  4.1.1  to  conclude  that  k  may  be  chosen  independently  of  X  We  shall  omit 
these  details  \~2 

In  general,  k a  could  be  too  large  to  allow  us  to  conclude  that  R *  is  a  k a-representation 
of  K*  However,  for  small  a  this  will  be  the  case  in  fact,  it  is  again  possible  to  find  a  <$i 
independent  of  A  such  that  R*  is  a  xa-representation  of  K*  whenever  a  <  5i 
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The  question  of  how  closely  the  model  K*  approximates  the  actual  intersection  of  K  with 
A  —  0  is  not  so  straightforward  as  the  first  question  because  we  have  not  established  a  stan¬ 
dard  for  compaung  polyhedra  with  different  face-topologies. 

One  final  question  concerning  accuracy  remains  to  be  addressed.  Namely,  given  an  a- 
representation  of  R  of  a  convex  polyhedron  K  and  linear  form  A,  how  do  we  estimate  suitable 
constants  5  and  k  'relative  to  the  given  A)  as  in  Theorems  4.4.1  and  4.4.2?  The  answer  is  to 
mimic  the  constructions  in  these  theorems  land  in  Theorem  4 .3.3  and  Appendix  C)  using  the 
data  of  R.  The  constructions  are  all  smooth  functions  of  the  numerical  data  so  that,  in  the 
end,  the  quantities  8  and  k  obtained  will  deviate  from  8  and  k  by  amounts  not  greater  than 
M  a,  where  iV/  is  a  constant  depending  only  on  the  face-topology  and,  perhaps,  some  gross 
boundedness  assumptions  As  long  as  8  >(M  +  l)a,  we  may  use  8'  —  M  a  to  estimate  8. 
Similarly,  we  may  use  a  +  M  a  to  estimate  a.  If  8'  S  (M  +  l)a,  then  we  must  reject  R  as 
being  insufficiently  precise. 

4.5  Iteration 

The  intersection  algorithm  of  §4.2  may  be  applied  to  the  output  representation  R*  and 
some  new  half-space  A*  ^  0  provided  that  R*  is  an  a*-representation  of  K*  far  suitably 
small  a*  We  have  seen  that  this  will  be  the  case  provided  we  have  been  able  to  choose  our 
initial  a  small  enough.  If  a  is  not  suitably  small,  R *  will  still  satisfy  conditions  (1)  -  (3)  of 
§4  1,  but  conditions  (4)  and  (5j  may  fail.  This  means,  essentially,  that  K*  has  certain  features 
that  are  too  small.  One  option,  then,  is  to  return  to  the  original  input  representation  R  and 
improve  its  accuracy  Another  option  might  be  to  single  out  the  undesirable  small  features  of 
K*  and  attempt  to  collapse  or  "shrink”  them  to  points,  much  as  we  collapsed  vertex  trees  in 
the  proof  of  Theorem  4.2.1  Such  a  shrinking  process  can,  in  fact,  be  carried  out.  but  it  has  the 
effect  of  moving  ail  the  vertices  of  the  polyhedron.  Thus,  once  some  features  are  ollapsed,  oth¬ 
ers  may  have  been  made  undesirably  small.  Moreover,  the  accuracy  of  the  modified  represen- 
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tation  may  not  be  as  good  as  that  of  R*. 


5.  Conclusions 

This  paper  explores  the  paradigm  of  structuring  code  so  that  all  conversions  from  numer¬ 
ical  data  to  symbolic  data  are  explicit  and  logically  independent.  The  paradigm  is  applied  suc¬ 
cessfully  to  the  problem  of  intersecting  convex  polyhedral  objects,  and  this  leads  to  an  algo¬ 
rithm  which  is  correct  in  the  sense  of  §1.3.  Since  correct  algorithms  are  clearly  robust  with 
respect  to  numerical-to-logical  conversions,  the  paper  produces  a  robust  algorithm  for  inter¬ 
secting  convex  polyhedra 

The  key  tool  used  in  the  paper  is  a  correspondence  between  stressed  planar  graphs  and 
convex  polyhedra.  which  allows  manipulation  of  convex  polyhedra  by  elementary  graph- 
theoretic  and  algebraic  means.  This  tool  is  used  to  deform  a  convex  polyhedron  so  that  certain 
of  its  vertices  may  be  placed  on  a  prescribed  plane.  It  may  also  used  to  remove  certain  small 
features  from  a  convex  polyhedron 

Prior  to  the  work  of  this  paper,  the  notion  of  attempting  to  make  conversions  to  symbolic 
data  independent  was  explored  in  the  implementation  of  an  algorithm  for  intersecting  noncon- 
vex  polyhedral  objects  [HHK2].  Although  the  algorithm  did  fail  on  certain  inputs,  its  numeri¬ 
cal  robustness  even  in  the  nonconvex  case  was  significantly  better  than  existing  algorithms. 
The  proof  techniques  of  this  paper  while  limited  to  convex  polyhedra  do  point  out  certain 
places  where  problems  are  likely  to  arise  in  intersecting  nonconvex  polyhedra.  For  example, 
the  underlying  gTaph  of  a  convex  polyhedron  is  always  3-connected,  and  essential  use  was 
made  of  this  property  in  various  proofs.  A  polyhedron  in  which  a  face  has  a  hole  (see  Fig.  5) 
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may  give  rise  to  a  disconnected  graph  Similarly  if  two  faces  intersect  in  a  disconnected  set 
(see  Fig.  6)  the  underlying  graph  will  not  be  3-connected.  This  suggests  that  when  either  of 
these  situations  occur  a  face  be  subdivided  so  that  the  graph  will  be  3-connected. 

The  topology  of  a  convex  polyhedron  places  only  minimal  geometrical  constraints  on  its 
vertices,  such  as  the  obvious  constraint  that  ail  vertices  on  a  face  lie  on  a  plane  In  the  case  of 
nonconvex  polyhedra,  however,  there  is  a  deep  interaction  between  geometry  and  topology.  In 
this  case,  for  example,  the  topology  can  force  three  vertices  to  lie  on  a  straight  line.  This 
allows  the  construction  of  a  topology  to  force  complex  geometrical  relations  between  vertices 
on  widely  separate  faces.  These  constraints  can  make  questions  concerning  the  incidence  of 
vertices  on  faces  logically  dependent,  and  thus  the  above  paradigm  fails.  However,  these  con¬ 
straints  seem  to  arise  only  when  the  intersection  of  two  faces  is  not  a  connected  region.  Thus 
checking  for  this  situation  and  partitioning  a  face  when  it  occurs  may  improve  robustness  sub¬ 
stantially 

Consider  two  rectangular  bars  crossing  (see  Fig.  7).  There  is  a  tendency  to  think  of 
ABCD  as  a  straight  line  that  is  the  intersection  of  Fi  and  F 2  However,  it  may  be  three  line 
segments.  In  fact  F2  and  F 2  may  not  lie  on  the  same  plane.  In  Figure  8,  F 1'  and  F 1  lie  on 
the  same  plane  By  making  use  of  this  construction  we  can  force  vertices  to  lie  on  a  straight 
line  and  hence  can  force  relations  between  various  vertices  by  means  of  theorems  such  as 
Pascal's  Theorem. 
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Figure  1:  Intersecting  two  polygons 


Figure  2:  Polyhedron  sitting  over  plane  and  the  graph 
of  its  projection 


Figure  3:  Examples  of  embedded  graphs,  with  and 
without  equilibrium  stresses. 


Force  vertices  on  different  faces  to  lie 
on  the  same  line 


Appendix  A:  Stressed  graphs  and  polyhedra 


This  appendix  provides  proofs  for  Theorems  3  3.1,  3.3.4,  and  3.6.1  (Corollary  A. 7,  Proposi¬ 
tion  A. 8,  and  Corollary  A. 10,  respectively).  We  begin  with  a  few  elementary  observations 
from  linear  algebra 

Every  u  e  R3  determines  a  linear  function  v*  :R3  — ►  R  by  the  rule  u*(w)  =  ww,  for 
every  w  e  R3  The  triple  product  [u,u,w]  =  det  (u,u,w)  may  then  be  expressed  in  terms  of 
the  cross  product  as  [u,v,iv]  =  ti*  X  w)*(u).  The  following  lemma  is  immediate: 

A.  1  Lemma:  \u  i ,  u2>  u3}i9  a  basis  of  R3  if  and  only  if 

/(u2  X  u3)*,(u3  X  X  u  2  )*} 

is  a  basis  of  (R3)  ;  where  (R3)  denotes  the  dual  of  R3,  i.e.,  the  vector  space  of  all  linear  func¬ 
tions  R3  -*  R.  Q 

Let  {U1.U2.U3}  be  a  basis  of  R3.  If  X  e  (R3)  ,  then  there  are  unique  constants 
c ! , c 2 ,  and  c  3  such  that 

X  =  C  [(U  2  X  U3)’  +  C2(U  3  X  U  1)'  +  C3(U[  X  u  2)'  (1) 

In  fact,  it  is  easy  to  see  that  c,  must  be  given  by 

c,  =  X(u,)/5  ,  1  =1,2,  and  3  (2) 

where  8  =  [u  t ,  u  2 ,  u3  ]. 

Consider  the  vector  space  A(H )  of  all  affine  functions  H  — »  R,  where  H  is  the  hyper- 
plane  of  R3  given  by  =  1. 
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A. 2  Lemma.  (R3)  is  isomorphic  to  A(H)  via  the  rule  A  —* ►  A  |  H  i  the  restriction  of  A  to  H). 

Proof:  The  mapping  A  — *  A|//  is  clearly  linear  Since  A|/f  =  0  implies  that  A  =  0  the 

maping  is  1-1.  It  is  onto  because  both  <  R3 )  and  A  kH''  *re  3-dimensional  [j 

Choose  three  noncollinear  points  P\,P2 •  and  P 3  in  W-  As  vectors  in  R3 ,  they  are 
linearly  independent.  Thus,  Lemmas  A.l  and  A, 2  yield  the  following  two  facts: 

A. 3  Corollary:  For  any  a  e  A(H),  there  are  unique  constants  c^^.and  c3  such  that 

a  =  c1(p2Xp3)'  +  C2(P3xPi)’  +  c3(pi  Xp2)' |H.  (3) 

In  fact,  setting  8  =  Ip  3 ,  p  2 ,  p  3  ],  it  is  easy  to  see  that  c,  must  be  given  by 

c,  =  a{pt)  8  ,  i  =  l,  2, and  3.  □  (4) 

A. 4  Corollary:  Suppose  a  and  a'  are  in  A(H).  and  a(pt)  =  a'(p.  for  1  =  1,  2.  Then 

there  is  a  unique  constant  c  such  that 

i 

a  -  a’  =  c(piXpi)'  \H.  □  (5) 

It  is  instructive  to  evaluate  c  in  this  case.  Choose  any  p.  e  H  not  collinear  with  P1P2 
Then,  (p  1  Xp2)  (p-)  =  [pi.p2,p*l  *  0.  Therefore,  after  evaluating  both  sides  of  (5)  at  p», 
we  may  solve  for  c: 

c  =  alp.)  -  a'lp*)  [P1.P2.P-I-  (6' 

It  follows  that  the  expression  on  the  right-hand  side  of  (6)  is  independent  of  p> . 

Now  recall  that  in  53  we  defined  A \PA(G)-*S  (G)  by  setting  A{ar}  =  where,  for 

each  edge  elJt  utJ  is  given  by 

Uij  =  e(r,  s)(as(p.)  -  ar(p-))  [pt , p , , p«  ].  (7) 

Here  fr  and  f\  are  the  faces  incident  to  elJt  and  e(r,s)  =  ±1  according  to  our  orientation 

convention,  as  described  in  $3.  The  preceding  paragraph  now  implies  the  following  result: 
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A. 5  Corollary  :  u>,,  is  independent  of  the  choice  of  reference  point  p*  i  1 


A. 6  Lemma:  For  any  i,  2  “i j  <P,  *P,>  =  0. 

;  =1 


Proof:  Let  f fo,  '  .  ft  denote  the  faces  of  Tip)  incident  to  p,,  in  counterclockwise  order 
around  p,.  It  will  be  convenient  to  write  r—*s  if  fs  is  an  immediate  successor  of  fr  in  this 
ordering.  Choose  any  r  =  1,2,  •  ■  •  ,  l,  and  s  such  that  r-»s.  Then,  choose  j  e  {1,2,  •  •  •  ,  n\ 
such  that  fr  n  fs  is  the  edge  connecting  pt  to  p^  Our  orientation  conventions  <*re  such  that,  if 
r—*s,  then 


.(P,i  =  i-l,  for  j>i, 
e(r's>~  1  l,  for  j  <i. 

If  J  >i,  then,  combining  (5),  < 6 > ,  and  i7),  we  see  that 


«,/P,  Xp7)' \H  =  -(a,  ~ar)  =  ar—as.  (8) 

If  j  <t,  then  j  and  i  in  formula  (7j  get  interchanged,  and  the  right-hand  denominator  changes 

sign,  as  does  e(r.s).  Thus,  we  get  equation  (8)  in  this  case  too. 

Now  sum  both  sides  of  (8)  as  r  =  l,2,  The  right-hand  side  telescopes  to  zero. 

Thus 


2  I#  =  °.  (9) 

j  =1 

from  which  the  desired  equation  follows  immediately  via  Lemma  A. 2.  □ 


A. 7  Corollary:  is  an  equilibrium  stress  on  Hp). 

Proof:  For  any  i  and  j,  p,  X  p;  =  (p,  —  p;)  X  p,.  Therefore,  by  Lemma  A  6, 

(  2  «.y(p,  -Pj))Xp,  =  2  «y(p,  Xp;)  =  0  , 

J  =  1  7  =  1 

a 

from  which  it  follows  that  2  ^i/Pi- Pj)  and  p,  are  linearly  dependent.  But,  the  former 

7=1 
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belongs  to  R2  =  (U,y,  0)  e  R3}  whereas  the  latter  belongs  to  H  =  {(x, v,  1)  e  R3}.  So  they 
can  be  linearly  dependent  if  and  only  if 

2  «,y(p.  ~Pj)  =  0  , 

'  =1 

as  required,  [j 

Theorem  3.3.4  asserts  that,  for  every  equilibrium  stress  there  is  a  P.A.  sequence 

{ar}  such  that  {u,^  =  A  (  ar\.  We  now  give  a  proof  of  this  fact.  The  proof  makes  use  of  some 
elementary  concepts  from  homology  theory.  For  a  cell  complex  K  and  real  vector  space  V,  we 
shall  assume  that  the  reader  is  familiar  with  the  concept  of  an  n-chain  on  K  with  coefficients 
in  V  and  the  corresponding  notions  of  boundary,  cycle,  and  homology.  We  suggest  [GH  ]  or 
[HY  ]  as  a  reference 

For  our  cell  complex  K,  we  use  the  2-dimensional  abstract  complex  given  by  Hp).  This 
consists  of  vertices  p(,  edges  etJ,  and  faces  fr,  with  the  incidence  relations  used  throughout  the 
paper.  For  our  vector  space  V,  we  use  A  {H),  the  vector  space  of  all  affine  functions  on  H. 

Typical  0-chains,  1-chains,  and  2-chains  will  be  denoted  by 

2  2  M«yl.  and  2  Crtfr], 

where  we  sum  over  all  vertices,  edges,  and  faces,  respectively,  and  at,blr  and  cr  are  ele¬ 
ments  of  A  [H). 

Given  any  equilibrium  stress  {u(;}  on  Hp),  we  now  define  a  specific  1 -chain  2 
by  setting 

btJ  =  -w./p.Xp/ \H. 

Note  that  btJ  =  —6,,,  for  all  i,j. 

A. 8  Lemma:  2  1  >8  a  1-cycle. 

Proof:  For  2  &iy[e,y]  to  be  a  1 -cycle  its  boundary  must  be  zero.  Formally,  the  boundary 

d(2  b,j[e,j])  is  a  0-chain  2a*tpt 3-  Thus  we  must  show  that  each  a,  =0.  Assuming  that  each 
edge  etJ  is  oriented  from  p,  to  p;,  the  usual  boundary  formula  yields,  for  i  =  1,2 . n, 
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a,  =  -  2  bu  +  2  bji  =  "  2  b:j  <  10> 

y  ='  =i  y=i 

=  2  u‘j(p<  xp.r  i#. 

j  =i 

Now  we  use  the  fact  that  {u,7}  is  an  equilibrium  stress: 

2  wiy(P<-Py>  =  0. 
y=l 

for  each  i  =  1,2,  ■  ,  n.  Therefore,  by  reversing  the  argument  of  Corollary  A. 7,  we  get 

2  Ui/p,  x-Pj)  =0, 
y=l 

which  implies,  using  (10)  above,  that  each  a,  =0,  as  desired.  Q 

So  far,  all  our  results  have  been  independent  of  the  specific  topology  of  K.  We  now  make 
use  of  the  fact  that  K  is  a  cellular  decomposition  of  the  2-sphere.  It  is  well  known  that  every 
1  -cycle  on  K  is  the  boundary  of  a  2-chain.  Thus,  by  Lemma  A. 9, 

2  b^e.j]  =#2  ar[fA),  (11) 

for  some  2-chain  2ar(/i-].  areA(H).  Here,  d  is  the  usual  boundary  operator. 

A.9  Proposition  For  every  equilibrium  stress  {w;y},  there  is  a  P.A.  sequence  {ar}  such  that 

=  A  {ar}. 

Proof:  Let  an  equilibrium  stress  be  given,  and  form  the  1-chain  aB  above. 

This  is  a  1-cycle,  so  as  in  (11),  it  is  the  boundary  of  a  2-chain  2 arlfrl-  We  ®how  that  {ar}  is  a 
P  A.  sequence  and  that  A{ar}  =  {u^}. 

Let  etJ  be  any  edge.  Formally,  the  coefficient  of  [ev]  in  d(2ar(/VJ)  »®  e(r ,  s)(ar  — as), 
where  fr  and  fs  are  the  faces  incident  to  etJ. 

Thus, 

«(r,s)(ar-as)  =  bt]  =  -  ul](pl  Xp;)' \H.  (12) 

Since  the  right-hand  side  vanishes  at  p,  and  p;,  ar  and  a„  coincide  on  etJ  Therefore,  { ar\  is  a 

P.A.  sequence 
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Furthermore,  if  both  sides  of  (12)  are  evaluated  at  p .  (which  is  noncoilinear  with  p,  and 
p,  i,  we  get 

uLJ  =  e(r ,  s)(as(p.)  -  ar{p.)).  [p,,py,p.  ], 

which  is  precisely  the  formula  used  to  define  A.  Therefore,  A{ar}  =  { utJ),  as  desired.  Q 

We  conclude  this  appendix  with  a  proof  of  the  convexity  criterion  (Proposition  3.6.1).  We 
begin  with  a  surface  S  sitting  on  the  plane  H  and  corresponding  P  A.  sequence  |ar},  with 
ao=0.  Set  }  =  A{ar}.  Clearly  a  necessary  condition  for  the  convexity  of  S  is  that  the 
plane  z  =  1  be  a  supporting  hyperplane.  Thus,  necessarily,  S  is  contained  either  in  the  half¬ 
space  z  S  1  or  in  z  S  1.  By  symmetry,  we  need  deal  only  with  the  former  case.  Recall  that, 
here,  we  always  assume  that  the  graph  T(p)  is  convex. 

A.10  Lemma.  Suppose  S  is  contained  in  the  half-space  z  S  l.  Then  the  following  are 
equivalent: 

(a)  S  is  strictly  convex. 

(b)  For  any  pair  of  adjacent  faces  fr  and  fs,  r  and  s  S  1,  and  for  any  qr  e  interior  (fr),  we 
have 

<M?r)  <  as(qr). 

Proof:  Let  fr  and  fs  be  adjacent  faces.  Select  qr  and  qs  in  the  interiors  of  frandfs  respec¬ 
tively.  Let  l(r,  s)  C  H  be  the  line  segment  connecting  qr  to  qs'.  that  is,  t(r,s)  = 
{tqr  +  (l-t)qs\0^t  SI}. 

(a)  (b).  Note  that  ar(qr)  *  as(qr),  otherwise  ar  =  as,  contradicting  strict  convexity. 

Thus,  if  (b)  fails  somewhere,  we  must  have 

ar(pr)  >  a^qr),  (13) 

for  some  r  and  s  as  in  (b).  Now  note  that  (13)  implies 

Sr(<?s)  <  as(qs).  (14) 

To  see  this,  restrict  both  ar  and  as  to  (( r,s ),  and  consider  the  function  values  of  ar(q)  and 
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as(q)  as  q  ranges  along  f(r,s)  from  qr  to  qs. 

It  follows  from  (13)  and  (141  that  every  value  of 

tar{qr)  +  (1 -t)as(qs),  0  <  t  <  1, 

is  strictly  greater  than  either  tar(qr)  +  (1  —  t)ar(qs)  =  ar(tqr  +  (1—  t)qs)  or 
ta^(qr)  +  (1  —  t)as(qs)  =  as(tqr  +  (1—  t)qs).  Therefore,  the  line  segment  in  R3  connecting 
( qr,ar{qr))t.  S  to  ( qs,a$(qs))eS  fails  to  lie  in  the  region  enclosed  by  S,  contradicting  convex¬ 
ity.  Thus,  (b)  must  hold. 

(b>  =*  (a).  The  idea  here  is  to  show  that  if  strict  convexity  is  violated,  there  exist  adja¬ 
cent  fr  and  fs  etc.,  for  which  (b)  fails.  To  find  such  fr  and  fs,  we  suppose  that  strict  convex¬ 
ity  is  violated  and  find  a  line  segment  l  connecting  points  on  S,  interior  to  distinct  faces,  such 
that  l  contains  a  third  point  on  or  above  S.  If  all  points  of  t  are  on  S,  then  S  has  coplanar 
adjacent  faces  which  can  serve  as  our  fr  and  fs.  Otherwise,  choose  qe(  whose  vertical  dis¬ 
tance  above  S  is  maximal  and  which  projects  vertically  to  an  interior  edge  in  G.  The  faces 
incident  to  this  edge  are  the  desired  fr  and  fs.  Q 

A.  11  Corollary:  S  is  strictlj:  convex  if  and  only  if 

utJ  <  0  ,  for  peripheral  edges  etJ, 

and 

u ij  >  0  ,  for  interior  edges  etJ. 

Proof:  This  follows  immediately  from  Lemma  A. 9,  together  with  formula  (7)  and  the  orienta¬ 
tion  conventions  described  in  S3.  O 

Remark:  A  similar  result  is  obtained  if  strict  convexity  is  replaced  by  convexity  and  strict 
inequalities  by  weak  inequalities. 
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Appendix  B:  A  Generalization  of  Theorem  4.3.3 

The  theorem  whose  proof  we  sketch  in  this  appendix  generalizes  Theorem  4.3.3  in  that  it 
applies  to  a  more  general  convex  polyhedron  K.  In  Theorem  4.3.3,  the  convex  polyhedron  K  is 
assumed  to  "sit  over  a  triangular  face"  in  the  plane  H,  as  described  in  §3.  In  the  more  general 
setting  of  this  appendix,  we  allow  K  to  sit  over  a  face  with  any  number  of  sides.  In  fact,  even 
more  generally,  we  consider  any  convex  polyhedron  K  with  no  faces  perpendicular  to  H,  and 
we  project  it  orthogonally  to  H.  Although  some  faces  of  K  may  have  overlapping  images 
under  this  projection,  each  individual  face  is  mapped  in  a  1-1  manner  to  H.  As  before,  the  1- 
skeleton  of  K  is  isomorphic  to  an  abstract  graph  T,  and  the  projection  of  this  1 -skeleton  into  H 
is  a  realization  of  T  denoted  Hp).  The  periphery  of  T(p)  is  a  convex  curve  in  H,  and  K  sits 
over  this  graph  as  before.  All  of  the  results  of  §3  relating  polyhedra  to  stressed  graphs  apply 
without  change  to  this  situation,  which  greatly  generalizes  that  considered  in  Theorem  4.3.3. 
For  example,  the  notion  of  a  p  a.  sequence  {ar}  is  defined  just  as  before,  as  is  the  notion  of  an 
equilibrium  stress  for  T(p).  And  the  correspondence  {ar}  <-->  {<jy}  is  defined  just  as  in 
the  more  restricted  setting  in  §3.  With  this  understood,  the  theorem  proved  in  this  appendix 
new  is  a  verbatim  restatement  of  Theorem  4.3.3. 

It  is  perhaps  worth  reminding  the  reader  at  this  point  why  such  a  generalization  is 
needed.  First  of  all,  we  start  with  a  plane  P  slicing  a  convex  polyhedron  K.  If  a  triangular 
face  of  K  has  vertices  a  reasonable  distance  from  P  on  each  side  of  P,  then  there  is  no  need  for 
this  generalization.  However,  this  may  not  be  the  case.  It  may  be  that  there  are  faces  with 
vertices  satisfying  the  condition  just  stated  but  none  that  is  triangular.  In  that  case,  we  nor¬ 
malize  just  as  before.  That  is,  we  choose  a  centra)  projection  from  a  point  on  P  near  the  said 
face  (and  outside  of  K),  and  project  onto  the  plane  of  this  face.  We  also  apply  a  projective 
transformation  so  that  this  plane  becomes  H,  P  becomes  the  plane  x  =  0,  and  central  projec¬ 
tion  becomes  orthogonal  projection  onto  H.  The  normalized  K  now  sits  over  this  non- 
triangular  face,  and  the  theorem  of  this  appendix  may  be  applied  to  it.  But  there  remains  a 
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further  possibility.  Namely,  each  face  of  K  may  be  effectively  on  one  side  or  the  other  of  P 
In  this  situation,  we  choose  an  edge  both  of  whose  endpoints  are  very  close  to  P,  and  we  per¬ 
turb  K  slightly  so  that  this  edge  is  moved  onto  P.  We  then  select  one  of  the  faces  containing 
this  edge  and  normalize  as  before.  In  this  case,  the  projection  point  should  be  chosen  on  P  and 
close  enough  to  the  face  so  that  the  resulting  peripheral  curve  consists  entirely  of  edges  of  the 
selected  face  and  those  of  its  neighbor.  Again,  the  theorem  of  this  appendix  may  now  be 
applied  to  the  normalized  K. 

We  now  start  with  K  normalized  as  just  described.  Let  p  denote  the  n-tuple  of  vertices 

of  the  projection  of  K  in  H,  and  let  X  and  Y  denote  the  corresponding  n-tuples  of  x-  and  y- 

coordinates,  respectively.  The  stress  constants  u:J  corresponding  to  K  determine  a  stress 
matrix  ft  as  in  §2,  i.e.,  a  symmetric  n  X  n  matrix  with  row-sums  0.  Let  e  and  f  denote  the 
number  of  edges  and  faces  of  K,  respectively.  Thus,  n  —  e  +  f  —  2.  The  spanning  tree  of 
rip)  has  n  1  edges.  Let  us  call  the  remaining  f  -  l  edges  of  Tip)  free.  If  we  assign  arbitrary 
stress  constants  to  these  edges,  and  we  choose  any  n-tuple  X  suitably  close  to  the  n-tuple  X, 
then,  as  in  Lemma  4.3.1,  we  may  solve  uniquely  for  the  remaining  stress  constants  so  that  x- 
equilibrium  is  satisfied.  That  is,  we  produce  an  n  X  n  symmetric  matrix  ft  satisfying 

a-x  =  o.  ns) 

Indeed,  since  ft  is  a  stress  matrix,  we  also  have 

ftl  =  0,  (16) 

where  I  =  (1,1,...,1)'  (cf.  the  proof  of  Lemma  2.3  1).  For  example,  if  X  =  X  and  the  "free" 
stress  constants  are  chosen  to  be  the  utJ  corresponding  to  free  edges,  then  we  obtain  the 
matrix  ft.  Let  s  denote  the  (/'•l)-tuple  of  these  free  stress  constants. 

Our  goal  is  to  find  an  n-tuple  Y,  depending  on  X  and  an  (/'-l>-tuple  of  free  stress  con¬ 
stants  s,  such  that,  writing  Y  =  y(X,s), 
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(a)  Y(X.s)  is  a  continuous  function  of  (X,s),  with  Y(X.s)  —  Y,  and 
bi  Q-Y  =  0.  (17) 

For  then,  by  restricting  close  enough  to  iX,s>,  we  can  force  Y  to  be  so  close  to  Y  that 

the  corresponding  graph  T(p)  is  convex.  Furthermore,  the  matrix  Q  can  be  forced  to  be  close 
enough  to  U  to  correspond  to  a  convex  stress  on  V(p),  which  by  ( 1 5)-<  1 7 )  is  an  equilibrium 
stress.  This  produces  the  desired  convex,  stressed  graph.  The  proof  now  concludes  exactly  as 
does  that  of  Theorem  4.3  .3. 

To  find  a  Y  satisfying  condition  i  1 )  above,  we  rely  on  the  following  result: 

B.l  Proposition:  rank  Q  =  n-  3. 

Proof  :  In  the  case  of  triangular  periphery,  this  fact  is  an  easy  consequence  of  the  results  in 
§2.  In  any  case,  since  Q  satisfies  U5M17)  for  X  =  X,  Y  =  Y,  it  is  clear  that  rank 
Q  S  n  —  3.  The  proof  of  equality  in  general  is  similar  to  the  proof  of  Theorem  5  of  [C]  and 
proceeds  by  contradiction.  If  rank  Q  <  n- 3,  then  there  exists  a  realization  Y(q)  of  T  in  R3 
projecting  orthogonally  onto  T(p),  having  3-dimeosional  convex  hull  Q,  and  having  {w,7}  as  an 
equilibrium  stress.  In  fact,  Q  is  the  convex  hull  of  the  vertices  q  i ,  tj 2 »  >4*  that  project  onto 

the  peripheral  vertices  Pi.P2>- ■■•Pk  of  Hp).  It  is  possible  to  find  two  such  vertices  that  are 
connected  by  an  edge  e  lying  in  the  boundary  of  Q  The  plane  through  e  perpendicular  to  H 
divides  Q  into  two  halves.  An  infinitesimal  rotation  of  one  of  these  halves  around  e  and 
towards  the  other  half  results  in  an  ifinitesimal  shortening  of  some  of  the  "interior"  edges  of 
T(q)  without  affecting  the  other  "interior"  edges.  This  shows  that  the  energy  function 
corresponding  to  T  (cf.  §2)  does  not  have  a  critical  point  at  q,  contradicting  Lemma  2.1.1.  I  I 

Now  we  relabel  vertices  if  necessary  so  that  the  n  —  3Xn  —  3  upper  left  corner  A  of  Q  is 
invertible  Then  the  same  will  hold  for  the  upper  left  corner  A  of  0,  provided  ( X,s )  is  suit¬ 
ably  close  to  ( X.s ).  Henceforth,  let  W0  and  Wj  denote  suitably  small  neighborhoods  of  X 
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and  s  in  R1  and  R'  1 ,  respectively,  and  restrict  (A'.s)  to  W0X  Then  we  may  write 


12  = 


A  B 
B‘  C 


A  (n  —  3i X(n  —3). invertible. 


A  simple  computation  shows  that 


nullspace  12  C  columnspace 


—  A  ~lB\ 
I 


where  1  is  the  3X3  identity  matrix.  Note  that  this  column  space  is  a  3-dimensional  subspace 
of  Ra  which  varies  smoothly  with  (X,s)  and  coincides  with  nullspace  12  when  (X,s)  =  (X,s). 
In  fact,  the  columnspace  =  nullspace  Q  if  and  only  if  rank!2  =  n  —3. 

Define  Y  =  Y(X.s)  to  be  the  orthogonal  projection  of  Y  into  the  columnspace.  Clearly 
this  satisfies  condition  (1).  Since  condition  (2)  is  equivalent  to  Y e  nullspace  12,  the  foregoing 
paragraph  shows  that  we  obtain  condition  (2)  whenever  rank.12  =  n  —  3.  Thus,  our  proof  will 
be  complete  as  soon  as  we  can  insure  this  equality 

To  formui<ue  tl  U  more  precisely,  it  is  convenient  to  recognize  that  the  spanning-tree  pro¬ 
cedure  of  Lemma  4  3  1  defines  a  smooth  mapping 

F 

WqXW!  —  S(n), 


where  S(n )  is  the  vector  space  of  all  nXn  stress  matrices.  The  dimension  of  S(n)  is 
—  n (n  —1).  Let  S(n,r )  denote  the  set  of  stress  matrices  of  rank  r.  It  is  easy  to  show  that 

mt 

S(n,r)  is  a  smooth  submanifold  of  S(n)  of  codimension  -^-(n  —  r)(n  —  r  —  1).  We  are 

interested  in  the  codimension-three  submanifold  S{n,n  —  3).  By  construction  and  Proposition 
B  1, 


F{X,s)  =  12  e  S(n,n-3). 


Let  M  denote  the  subset  F  1  (S  in,  n  —  3))  C  Wq  X  W  j ,  and  let  ir:  VV0  X  W  t  — *W 0  be  the  pro- 
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jection  maf  it\X,s)  =  X  .  The  following  is  our  main  technical  result. 

B.2  Proposition:  Assume  that  the  graph  Tip)  has  a  triangle  li  e..  a  simple,  closed  edge  path 
consisting  of  three  edges).  Then 

(a)  There  exists  a  neighborhood  V  C  WqX  W[  of  (X,s),  such  that  M  n  Visa  smooth 
manifold  of  dimension  n  +  f—  2. 

(bi  There  exist  open-ball  neighborhoods  Bq  around  X  in  W0  and  B  x  around  0  in  R'  ~2, 
together  with  a  diffeomorphism  h:B0XB  i  — »M  n  V  satisfying 

■nhixi,  ■  ■  ■  ,xn,z\,  ■  •  •  ,zf  _2)  =  (xj,  ■  ■  ■  ,xn)  , 

for  all  (x ! ,  •  •  •  ,xn,zi,  ■  ■  , zt-  -  2)  in  fio  1  • 

Let  8  be  the  radius  of  B  o,  and  suppose  \X  —  X]  <  8.  Define  s  by  the  equation  (X,s)  = 
h(X,  0).  Then  s  depends  smoothly  on  X,  and  F(X,S )  has  rank  n-3  by  construction.  Thus, 
assuming  Proposition  B.2  and  the  triangle  hypothesis,  we  are  done. 

We  now  indicate  the  proof  of  Proposition  B.2  and  then  how  to  get  rid  of  the  triangle 
hypothesis.  Basically,  it's  a  simple  application  of  the  Inverse  (or  Implicit)  Function  Theorem. 
The  hard  part  is  computing  that  a  certain  derivative  is  non-zero.  We  shall  give  a  few  details. 

Rewrite  Q  as 

A  B 1  B 2  A  (n  —3) X(/i  —3),  invertible 

B[  Cn  Cl2  B  \  {n  —  3)X2 

B‘2  C‘12  C 22  B2  U-3)X1 

B.3  Lemma:  rankQ  =  (n  — 3)  +  rank(C u  —  B[A  ~lB  1).  □ 

Thus,  (X,s)  e  M  provided  rank(Cn -B\A  ~lB \)  0.  Now  Cu-B\A~lB  1  is  a  2  X  2 

symmetric  matrix  of  rank  £  1  (because  1  and  X  belong  to  nullspace  Q)  Such  a  matrix  is  zero 
«=>  it  has  zero  trace.  So  define 
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g{X,s)  =  trace{Cu  -B[A  ~lBl). 

This  gives  a  real-valued  function,  g  :Wo  X  W  i  — »R.  such  that 

F  ~l(Sin,n  -3))  =  M  =  g  _1(0). 

The  advantage  of  g  over  F  is  that  g  is  real-valued,  and  the  kind  of  regularity  result  that  we 
want  is  easier  to  obtain  for  g. 

Now  for  any  XeWo,  define  gx  to  be  g  restricted  to  the  "slice"  {A]  X  W  x  .  more  precisely 
gx'-W  i— »R  is  given  by  g\(s)  =g(X, s).  We  then  have  the  following  result: 

B.4  Proposition:  Assume  that  the  graph  Hp)  has  a  triangle.  Then  s  is  a  regular  point  of  gx 
That  is,  d(gx)-;  0.  Consequently,  (X,s)  is  a  regular  point  of  ir  |  M  :  M— ^Wq. 

Proposition  B  2  follows  easily  from  this. 

The  second  assertion  of  Proposition  B.4  follows  from  the  first  by  standard  arguments  of 
differential  topology.  Therefore,  it  remains  to  prove  d(gx  )j  *0.  To  show  this,  we  need  only 
compute  this  derivative  in  one  judiciously  chosen  tangent  direction.  Such  a  tangent  direction 
is  an  appropriate  ( f  —  l)-tuple  of  free  stress  constants.  Now,  with  our  triangle  hypothesis,  we 
are  able  to  choose  such  constants  so  that  the  corresponding  (infinitesimal)  stress  matrix  has  at 
most  three  non-zero  rows  and  three  non-zero  columns.  Thi*  greatly  simplifies  the  computation 
and  is  the  only  reason  for  the  hypothesis.  The  computation  proceeds  by  cases  according  to  the 
distribution  of  the  selected  rows  and  columns  in  the  matrix  There  are  seven  cases.  We  shall 
go  through  one  case  here  for  illustration. 

We  make  a  preliminary  observation.  A  stress  on  T  satisfying  ^-equilibrium  with  respect 
to  A'  is  uniquely  determined  by  the  stress  assignments  to  the  free  edges.  Thus,  instead  of 
selecting  an  assignment  to  the  free  edges,  we  select  an  entire  stress  This  is  notationally  con¬ 
venient,  since  then  we  do  not  have  to  distinguish  between  stress  constants  that  are  free  and 


those  that  are  not  free 
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In  the  case  that  we  consider  here,  the  triangle  of  the  hypothesis  has  do  vertical  edges 
and  has  vertices  p„  -2-Pn  -i-Pn  We  choose  a  stress  satisfying  x-equi librium  as  follows.  To 
*11  edges  not  in  the  triangle,  assign  stress  constant  0.  The  edges  pn  -iPn  -i-Pn  -iPn-Pn  -2 Pn 
get  the  stress  constants  1,  a,  b.  respectively,  where  a  and  b  are  determined  by  the  x- 
equilibrium  requirement.  It  is  not  hard  to  check  that  (1  +a)(l  +  b)  =  1,  again  because  of  x- 
equilibrium.  In  this  stress-direction,  one  calculates  that  the  differential  d(gx h  has  the  value 

trace  [  1 +/>  lfl  ]  =  a  +  b  +  2, 

which,  by  the  preceding  relation  on  a  and  6,  has  absolute  value  greater  than  two.  Thus,  Pro¬ 
position  B.4  is  verified  in  this  case.  Complete  details  will  appear  in  [K2] 

This  completes  our  discussion  when  T(p)  has  a  triangle. 

In  case  T(p)  does  not  have  a  triangle,  it  must  have  an  order-three  vertex  (in  fact,  at  least 
eight  of  these).  Perform  a  "slice"  at  this  vertex,  replacing  it  by  a  small  triangle.  Carry  out 
the  argument  for  the  modified  graph.  Then  restore  the  vertex  by  reversing  the  slice  pro¬ 
cedure. 


